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A novel general formalism for the maximal symetrization and reduction of fields (MSRF) is pro- 
posed and applied to wavefunctions in solid state nanostructures. Its primary target is to provide 
an essential tool for the study and analysis of the electronic and optical properties of semiconductor 
quantum heterostructures with relatively high point-group symmetry, and studied with the k ■ p 
formalism. Nevertheless the approach is valid in a much larger framework than k ■ p theory, it is 
applicable to arbitrary systems of coupled partial differential equations (e.g. strain equations or 
Maxwell equations). This general MSRF formalism makes extensive use of group theory at all le- 
vels of analysis. For spinless problems (scalar equations), one can use a systematic Spatial Domain 
Reduction (SDR) technique which allows, for every irreducible representation, to reduce the set of 
equations on a minimal domain with automatic incorporation of the boundary conditions at the 
border, which are shown to be non-trivial in general. For a vectorial or spinorial set of functions, the 
SDR technique must be completed by the use of an optimal basis in vectorial or spinorial space (in 
a crystal we call it the optimal Bloch function basis (OBB)). The full MSR formalism thus consists 
of three steps : 1) explicitly separate spatial (or Fourier space) and vectorial (spinorial) part of the 
operators and eigenstates, 2) choose, according to the symmetry and well defined prescriptions (e.g. 
specific transformation properties), optimal fully symmetrized basis for both spatial and vector (or 
spin) space, and 3) finally apply the SDR to every individual scalar ultimate component function. We 
show that with such a formalism the coupling between different vectorial (spinorial) components by 
symmetry operations becomes minimized and every ultimately reduced envelope function (UREF) 
acquires a well-defined specific symmetry. The advantages are numerous : sharper insights on the 
symmetry properties of every eigenstate, minimal coupling schemes (analytically and computatio- 
nally exploitable at the component function level), minimal computing domains. The formalism can 
be applied also as a postprocessing operation, offering all subsequent analytical and computationnal 
advantages of symmetrization. The specific case of a quantum wire (QWRs) with Cz v point group 
symmetry is used as a concrete illustration of the application of MSRF. 

PACS numbers: 73.21.-b, 78.67.-n 



I. INTRODUCTION 

The study of low dimensional solid state nanostruc- 
tures is a very interesting and promising domain. Indeed 
a good knowledge of electronic and optical properties 
of nanostructures like metallic [l| or semiconductor [2[ 
nanostructures, or photonic crystals [|[, is essential for 
many applications in advanced lasers, photonics and te- 
lecommunications. New truly quantum applications like 
quantum cryptography also make extensive use of semi- 
conductor quantum heterostructures like quantum wells, 
quantum wires and quantum dots [H [j| . The quality of 
semiconductor quantum wires and dots have been exten- 
sively improved during the last fifteen years, and one is 
now able to produce high quality structures with higher 
and higher point-group symmetries (e.g. Cq v quantum 
dots [a, [7[). In such a case a group-theoretical approach 
is usually the most powerful tool for describing the ef- 
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fects resulting from symmetry on the electronic states, 
their optical properties in particular. However the pro- 
blem is rather complicated since in heterostructures one 
must take into account both the underlying microscopical 
crystalline structure and the mesoscopic heterostructure 
confinement potential. 

The theoretical study of low symmetry effects in semi- 
conductor heterostructures (like quantum wires with C s 
symmetry, e.g. T- and V-shaped quantum wires [H, Q) 
is already well developed and has lead to funda- 
mental conclusions regarding their electronic and opti- 
cal properties. First electronic and excitonic states can 
be labelled with respect to their characteristic transfor- 
mation properties under symmetry operations, second ri- 
gorous and important selection rules were readily obtai- 
ned on the basis of such a classification, useful even in a 
low symmetry case [T3, El ■ The effects of lateral confi- 
nement to the polarization anisotropy were largely stu- 
died p^4l6j . However, it should be pointed out that up 
to now only very little work has been devoted to higher 
symmetries (called here HSH : High Symm etry Hetero- 
structure), for example C$ v structures [It], Il8j . or even 
Cq v [lj| . A particularity of HSH is to allow the existence 
of symmetry-induced degenerate eigenstates, related to 
irreducible representations (irreps) with dimension grea- 
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ter than one [2JJ,|2l|, and which display a much more com- 
plex behavior under symmetry operations. Symmetry- 
induced degenerate eigenstates play an important role in 
the generation of entangled photon pairs from quantum 
dots [5|. 

The electronic structure of semiconductor heterostruc- 
tures are very often studied in the frame of the k ■ p en- 
velope function approach for heterostructures 22], with 
at least four bands when describing the valence band. 
In such a frame the different envelope functions (com- 
ponents of the spinorial eigenstates) become entangled 
under symmetry operations : their shapes are therefore 
mutually coupled, and their behavior is complex. Up to 
now there has been very few attempts to try to use an 
Optimal Bloch function Basis (OBB). In [H| we tried 
to rely on the concept of an 'Optimal Quantization Axis 
(OQA) direction' (A Bloch function basis which diagona- 
lizes the component of angular momentum in the chosen 
optimal direction). In fact we shall show in the following 
that such a method is of limited interest : it is optimally 
adapted only in very low symmetry cases like C s struc- 
tures ! For Quantum Wires (QWR) with a higher sym- 
metry group, the previously defined OQA direction may 
only be an improved choice, not the optimal choice. 

In this paper, we propose the optimal and systematic 
solution to this problem : a general Maximal Symmetriza- 
tion and Reduction (MSRF) formalism, perfectly adap- 
ted to the study of scalar or spinorial HSH problems. We 
will show how to find true OBBs which minimize the 
coupling between envelope functions, and which truly 
maximize their individual symmetry. Moreover we will 
show how to systematically compute the whole solution 
on a reduced minimal domain. With the improvement of 
growth techniques increasingly high symmetries can in- 
deed be produced enhancing the need for novel tools al- 
lowing to fully take into account the symmetry properties 
and to significantly simplify, theoretically and numeri- 
cally, the understanding and the computation of electro- 
nic and optical properties. We originally developed the 
MSRF formalism to study a Cs v quantum wire (QWR) , 
and therefore we shall often use it as an example, but 
one should stress that the method is general, applicable 
to other groups and to widely different cases since in 
fact its possible scope of application is much wider : it 
could be applied in its full generality to arbitrary ten- 
sorial fields obeying a set partial differential equations 
characterized by any given point group symmetry. In 
particular the method is independent of the number of 
coupled bands kept in the problem, and independent of 
the specific terms kept, provided the global symmetry is 
conserved. For example we could easily take into account 
interface terms or strain terms in an eight-band 

approach 25, 26j , but in case of strain we would also need 
to treat in the same way the elasticity equations for the 
strain tensor (this could be done most conveniently by 
post-symmetrization of the elasticity calculation, see end 
of Section l*VT)) . The MSRF method is also independent 
of dimension, it applies equally well for two or three spa- 



tial dimensions, i.e. to QWR or to quantum dot (QD) 
heterostructures with a given point-group symmetry. 

Let us now shortly present the heart of the MSRF 
formalism, which is threefold. First for every quantum 
states one performs an explicit separation of the spatial 
character (3D orbital motion, eventually treated in Fou- 
rier space) and of the field character (e.g. spinorial). Se- 
cond one selects optimal fully symmetrized bases, both 
for spinorial space (the Optimal Bloch function Basis 
(OBB)) and orbital space, which minimize the coupling 
between different spinorial components. These two ingre- 
dients allow to obtain every spinorial component of the 
field as a sum of symmetrized scalar functions of spa- 
tial coordinates. Third for every irrep one identifies mi- 
nimum sets of independent parameters (orbital reduced 
domains) which form the systematic Spatial (or Fourier) 
Domain Reduction technique (SDR) and which allow to 
obtain reduced Hamiltonians with respect to reduced do- 
mains, and systematically minimally reduce the compu- 
ting requirements. 

The advantages of the new MSRF formalism are mani- 
fold. Indeed besides the possibility of performing SDR we 
shall show that there are also many advantages from the 
analytical point of view : first the Hamiltonian operator 
usually takes a simpler form in the adapted fully sym- 
metrized basis (OBB), second the spinorial components 
of eigenstates (as well as the components of any operator 
in the spinorial basis) can be treated in a similar way 
and easily decomposed into fundamental parts to which 
single group irreps can be associated (and for which " sub- 
selection rules" can be applied at an intermediate calcu- 
lational level). In this way particularly simple analytical 
expressions can be obtained for certain operator matrix 
elements, which allow to find, for example, new analytical 
ratios in the polarization anisotropy that were previously 
unnoticed in the numerics [27]. Further insight can also 
be gained from the fact that this symmetry-based tech- 
nique simplifies the expression of coupling matrix ele- 
ments. Most notably weak symmetry breaking mecha- 
nisms can be understood more deeply at the analytical 
level [28[ . From the numerical point of view, the systema- 
tic SDR technique will enable one to solve independently 
for every irrep on a minimally reduced solution domain. 
The SDR technique not only allows to find the boundary 
conditions at the boundary of such a domain, shown be- 
low to be non-trivial, but it also allows to eliminate the 
need to explicitely care for them ! It should be pointed 
out that the MSRF method can also be used as a post 
symmetrization technique on numerical results obtained 
without taking into account any symmetry. In such a case 
it not only allows to classify all eigenfunctions and sym- 
metrize them within the OBB, and benefit of an in-depth 
symmetry analysis, but in all subsequent computations 
symmetrized wavefunctions on reduced domain can then 
be used, which may still represent a further significant 
potential gain. 

Let us now detail the necessary procedures of the pro- 
posed MSRF technique. 
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For scalar problems, like the single band k ■ p spinless 
conduction band Hamiltonian, it reduces to our syste- 
matic SDR procedure. The SDR procedure involves two 
fundamental steps : first the spatial domain must be de- 
composed into a minimal number of disjoint sub-domains 
which map into each other through symmetry operations 
(including borders as separate domains), second the set 
of domains and wavefunctions must be projected on the 
relevant irreps. This last step allows to find the critical 
geometrical features of all states by identifying for every 
irrep the minimal independent parts of any function of a 
given symmetry, to which a corresponding reduced sub- 
domain can therefore be associated. At this stage non- 
trivial boundary conditions can be derived if necessary. 
The same procedure can then be carried out for all rele- 
vant functional operators like the Hamiltonian. The re- 
duced Hamiltonian reflects directly the coupling between 
different sub-domains and does incorporate automatically 
the restrictions implied by non-trivial geometrical boun- 
dary conditions. 

For spin dependent problems, like the 4 x 4 k ■ p Luttin- 
ger Hamiltonian describing the valence band in diamond 
semiconductors, or the much used eight-band k-p Hamil- 
tonian [26[, the OBB basis functions must be first found. 
They transform like an irrep of the (double) group, and 
allow to block-diagonalize the corresponding matrix re- 
presentation of the double-group. By choosing the OBB, 
one enforces a minimal coupling between different spi- 
norial function components under symmetry operations. 
Since every component can then be decomposed in a 
simple way into scalar envelope functions labelled with 
single groups irreps, the SDR technique is applicable 
to every component, and reduced Hamiltonians can be 
found for every double-group irrep. 

In section |TT] we shall recall in more details essential 
results obtained in low symmetry heterostructures which 
are needed to understand the HSH challenge and establi- 
shing the basis for the development of the MSRF forma- 
lism. In Section Hill transformation laws in both ordinary 
space and spin space are studied, together with funda- 
mental group-theoretical results needed for the definition 
of the fully symmetrized OBB basis and for the separa- 
tion of spinorial and spatial parts. These goals are attai- 
ned in Section HV|) . In the next two Sections, we first de- 
velop the SDR formalism for single group spinless (scalar) 
functions (Sec.|V|, and in the second we apply the SDR 
to the symmetrized envelope functions created by the 
OBB (Sec. IVip . In addition we show how the technique 
leads to reduced Hamiltonians, and how it can be used as 
a post-symmetrization technique. Finally, in Section [VIII 
we demonstrate how selection rules can be applied at an 
intermediate level - a specific feature of the MSRF for- 
malism - to compute the matrix elements of operators. 
As a result we find novel strong analytical results (pola- 
rization anisotropy in C^ v structures) which can be in- 
terpreted with the help of the generalized Wigner-Eckart 
theorem for point groups. Finally, in Section rVIIIi an out- 
look is given on other symmetry groups (the hexagonal 



Cq v group, an approximate zone-center symmetry group 
D 3 h , the (commutative) C n , n £ N subgroup of the rota- 
tion group and the C s group) . The potential of the MSRF 
formalism for different problems is shortly described in 
Section HXl as well as its relationship with the most close 
works found in the litterature on heterostructures which 
exploit symmetry. 



II. ENVELOPE-FUNCTION THEORY OF LOW 
SYMMETRY HETEROSTRUCTURES AND THE 
HSH CHALLENGE 

The MSRF formalism developed in this paper is built 
upon specific techniques previously developed for low 
symmetry heterostructures. The main goal of this sec- 
tion is not only to introduce the basic envelope function 
k ■ p Hamiltonians for the conduction and valence band 
that will be used throughout the paper, but also to re- 
call fundamental results on low symmetry heterostruc- 
tures [l(| which are at the origin of MSRF. We will also 
show explicitly the limitations of these techniques which 
do not allow to reach maximal symmetrization for higher 
symmetry heterostructures, which thus represent a main 
challenge. These results will form an essential basis for 
the systematic study of transformation laws developed in 
Section iHll and the development of the cornerstone of the 
MSRF in Section [IVJ 



A. Introduction to envelope function models 

Multiband k ■ p Hamiltonians are extensively used for 
the study of the electronic structure in semiconductor 
heterostructures [22| . Such models allow to introduce all 
the relevant physics close to a high symmetry point of 
the band-structure, whilst keeping maximum simplicity. 
For many applications in III-V zincblcndc semiconduc- 
tors like AlGaAs/GaAs it is possible to treat separately 
the conduction and the valence band problems. Never- 
theless the method presented in the following is genera- 
lizable to more complex multiband schemes which treat 
simultaneously the coupling between these bands, like 
eight-band or fourteen-band k ■ p Hamiltonians. 

Let us also assume a heterostructure translation with 
translation invariance in some spatial directions and 
which are confining in the remaining directions, like 
quantum wells or quantum wires (see Fig. [TJ, and split 
the coordinate system into ry and rj_ respectively. 

For the isolated conduction band the k ■ p approach 
gives rise to the simple effective-mass model when one 
ignores spin-splitting terms. In the case of QWRs or QWs 
the conduction band Hamiltonian operator H is defined 
by its action H [tp] on any electron state ip. In the r repre- 
sentation, this amounts to apply the following differential 
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operator 

ti 2 1 
2 m(r_i_) 



2 m (rj_) 



+ V c (r x ) 



on the electron wavefunction 7/>(r_i_)) which is a simple 
scalar function of position. Here the perpendicular gra- 
dient Vi is of course related to the confined directions, 
and translation invariance of H implies a translation 
invariant confining potential V c (rj_) and effective mass 
m(rj_), both independent of ry, and the appearance of 
the corresponding good quantum number ky , which can 
be interpreted as the electron momentum in the free 
directions. Note that we have used for convenience in 
Eq. (U) the somewhat clumsy notation H(r±,'k\\), which 
comprises differential operators like Vj_ and variables like 
rj_ and k\\ , keeping in mind that, in the following, even- 
tual transformations on the argument rj_ of iJ(rj_,ky) 
must also be applied consistently to Vj_. Let us also 
denote ifa., (rj_), the eigenstates of if(rj_,kn) associated 
with the eigenvalue E^.. . 

For the valence-band let us use the minimal four band 
Luttinger Hamiltonian [2^, [29l - [3l| which is required when 
one wants to have a good estimate of the optical polari- 
zation anisotropy. Such model also provides a fairly good 
description of the QWR valence subband energy disper- 
sion close to the T zone-center. 

In its original form the bulk Luttinger Hamiltonian can 
be written as (2!| [3(| 

Hl = -— {|(7i + |72)fc 2 -7 2 E fc ^ 
mo 12 2 



le/]kikj(JiJj + JjJi) 



(2) 



where O denotes all cyclic permutations of the indices 
(i j € {1,2,3}) of the 4x4 matrix representations 
of the components of quasi-angular momentum Ji , i = 
1, 2, 3. The corresponding envelope function Hamiltonian 
for the QWR valence band reads : 



H(T±,k)=H L (T±,k)+hV v (T±) 



(3) 



where l4,V v (r±) is the diagonal part incorporating the 
effect of the heterostructure confinement potential, here 
k = k|| is the wave vector along the wire, and H L (rj_, k) 
is the kinetic part as given by the bulk Luttinger Hamil- 
tonian where rj_ — > — iVj_. As a result the kinetic part 
can always be put in the following standard form |3l| : 



also involve spatially-dependent Luttinger parameters 
7i( r J-),* = 1,2,3 corresponding to the spatial composi- 
tion dependance of the heterostructure. The eigenstates 
of H (rj_, k) can be considered as 4D spinorial fields (re- 
lated to the j = 3/2 Bloch function basis at the top of 
the valence-band), and each component of this field is 
a /c-dependent scalar field of rj_, the so-called envelope 
functions. 

At this stage a few fundamental comments are in or- 
der. First the actual form of the p, q, r, s coefficients is 
a function of the Bloch function basis chosen. The stan- 
dard form (3l| . corresponding to the main crystal direc- 
tions [100], [010], [001], is obtained using standard matrix 
representations [32] of Ji,i = 1,2,3. However for hete- 
rostructures whose main symmetry elements differ, one 
should usually define a rotated cartesian frame with a 
new z-axis oriented along adapted directions, e.g. see [3lJ 
for [hhk] directions, such that the corresponding implicit 
Bloch function basis diagonalizes the new J z component 
of the quasi-angular momentum. 

It is however important to point out that the shape of 
each envelope function is basis- dependent, even though 
the expectation value of any physical quantity remains 
basis-independent ! This point, which will be illustrated 
in the next subsection, was not explicitely recognized for 
a long time, and complicates a lot the question of the 
symmetry of the envelope functions in the presence of 
valence band mixing. 



B. The effects of symmetry on the envelope 
functions 

It is well known [33| that the existence of a symme- 
try group for an Hamiltonian allows to classify each kind 
of eigenstates, to deduce their degeneracies, and to spe- 
cify their possible transformation properties, however in 
the case of a spinorial set of envelope functions, the 
analysis of the symmetry of every individual envelope 
and its transformation properties has not been addres- 
sed and requires much more work (Section IIII[) . Let us 
give here an intuitive justification for this extensive ef- 
fort. We shall first discuss the effects of a single symmetry 
plane, already quite well-known (see [10( , and subsequent 
works [ill [34Tj3a |V Then we shall enlighten in detail the 
difficulties involved in applying these results, based on ra- 
ther elementary concepts, to higher symmetries. The dis- 
cussion will both develop intuition and a new understan- 
ding of the difficulties involved, demonstrating the need 
for MSRF. 



Hr 



p + q —s r 
— s + p — q r 
nir, I r+ p — q s 

r + s + p + q 



(4) 



In this equation the p, q, r, s matrix elements are fc- 
dependent partial differential operators obtained from 
the Luttinger quadratic polynomials ([2"]). Eq. (jj) may 



1. A low symmetry case : C s symmetry 

Let us now assume a heterostructure with the sim- 
plest symmetry, like the typical V-shaped QWR shown in 
Fig. [TJ with a single symmetry plane. The corresponding 
symmetry point group is C s . Let us take the coordinate 
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system such that rj_ = (y, z). The symmetry of the hete- 
rostructure with respect to the plane implies that V c (rj_) 
and m(rj_) are both invariant with respect to the sym- 
metry plane operation a : rj_ i— > cr(rj_) = (y,—z). The 



the cr-symmetry plane 



+< ) (r ± ) 
-^(rx) 



(5) 



A y [ooi] 




Figure 1: View in perspective of a typical V-shaped 
Al x Ga\- x As QWR, with coordinate system. The central part 
is typically pure GaAs, whilst the surrounding bulk part is ty- 
pically Al x Ga\- x As with x ~ 30%. 



wavefunction profiles of the first two eigenstates of the 
stationary Schrodinger equation are shown on Fig. [2] It is 
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Figure 2: Contour plots of conduction band envelope func- 
tions of the first two electronic states in a typical V-shaped 
QWR 

(a) Even function (ground state) ; (b) Odd function (first ex- 
cited state). 

The V-shaped QWR GaAs potential well is shown with a 
doted line, the vertical quantum well with x m 20% is also 
visible. 



easy to show that the symmetry of the structure implies 
that the eigenstate symmetry can be labelled as being 
either "even" or "odd" . Indeed we see on Fig. [5] that 
the ground state wavefunction ip^i r ±) is strictly even 
with respect to the symmetry plane a, whilst the first 
excited state wavefunction tp^\r±) is strictly odd. The 

higher states will all display one character or the other 
unless there would be an accidental degeneracy, which 
would then allow an accidental mixing within the dege- 
nerate subspace. To summarize : the eigenfunctions all 
obey one of the following transformation rules under 
the operation of reversing a wavefunction with respect to 



It is also obvious that relations (|5|) do translate into strin- 
gent conditions for the properties of the wavefunctions on 
the symmetry axis : 



£>W£,W = 0) = 



0) = 



(6) 



Such relations are very useful because they can be used as 
boundary conditions to reduce the domain of solution on 
the left or right halfplane, which is therefore the natural 
reduced domain of solution of the stationary Schrodin- 
ger equation in this case. In most solution schemes, e.g. 
real space methods like finite element (FE) or finite diffe- 
rences (FD) approaches, it is easy to obtain odd and even 
solutions separately by solving two times the eigenpro- 
blem with different Dirichlet /Neumann boundary condi- 
tions on the symmetry axis boundary. Let us shortly have 
a group-theoretical approach : the even and odd wave- 
functions shown in Fig. [5] correspond respectively to the 
A' and A" irreducible representations (irreps) of the C s 
group. A' and A" are simply new group-theoretical labels 
meaning "even" or "odd" , and no further insight arises. 
Nevertheless our considerations related to the symmetry 
of conduction band wavefunctions in C s symmetry - des- 
pite their trivial aspect - will prove important to better 
appreciate the differences occurring for valence-band en- 
velope functions. 

The case of the valence band eigenstates in C s sym- 
metry is much more complex because of their (four- 
dimensional) spinorial character, this is why point group 
theory will immediately become an unvaluable asset. It 
tells us [2l| that the spinorial eigenstates bare two la- 
bels again, but this time corresponding to the double 
group irreps l Ei/ 2l i — 1,2 instead of A' and A". The 
spinors bearing these labels display a much more subtle 
and disturbing behaviour : there is not any more any 
simple intrinsic symmetry for each of the envelope func- 
tion components ipk,m(y,z) [where fc = |kn |]. Their indi- 
vidual symmetry may indeed depend on the basis chosen, 
and, although there are two kinds of eigenstates, there is 
no such simple geometrical interpretation of the symme- 
try of such states by words of every-day life like "even" 
or "odd" . The role of their more complex label l E x / 2 is 
precisely to convey the nature of these more complex 
transformation laws under mirror symmetry. 

Let us now illustrate this behaviour with the valence 
band eigenstates of the V-shaped QWR shown in Fig. [T] 
With the standard Bloch function basis used in the early 
works on the subject [H, HH, the ground state envelope 
functions are those of Fig. [3J where one clearly sees that 
none of the envelope functions is either perfectly sym- 
metric or antisymmetric with respect to the symmetry 
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plane ! This "standard" choice of Bloch function basis 
was at the time guided by the fact that the shape of 
a V-shaped QWR is close to a deformed quantum well, 
and it was indeed reasonable since many qualitative fea- 
tures of the optical absorption spectrum could be un- 
derstood [H, IHI on the basis of "quantum well light and 
heavy holes". This is why this basis is the one diago- 
nalizing the component of pseudo-angular momentum J 
aligned with the [001] crystalline direction (note : with 
the choice of labels as in Fig. [T] this actually correspond 
the vertical direction y, i.e. J y ). 
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Figure 3: Contour plots of valence band envelope functions 

tp k ^ 2 {y,z) of the ground state spinor (symmetry 1 -E 1 / 2 ) in a 
typical V-shaped QWR, when the Bloch function basis diago- 
nalizes J y corresponding to the [001] crystal direction (same 
QWR as in Fig.HJ). (a-d) m = 3/2, 1/2, -1/2, -3/2 function 
components respectively {Jy). 

Even if the envelope functions of Fig. [3] are not sym- 
metric, intuitively one still expects some symmetries in- 
duced by the QWR symmetry. Indeed a closer analytical 
look reveals that there are still some symmetry relations 
linked with 1 2?i/2 or 2 Ei/ 2 eigenstate. They can be for- 
mulated as follows : 



•^jk,-m(V»-») 



(7) 



Clearly such symmetry relations, which only couple ±m 
envelope functions, cannot enforce the individual sym- 
metry of every envelope functions in the spinor, and are 
nevertheless awkward from the numerical point of view 
since they do not allow to reduce the domain of solution 
on the half-plane as in the spinless case ! 

The clue to this problem was found in [10] by choo- 
sing a different Bloch function basis which diagonalizes 
the component J z oriented along the [110] crystalline di- 
rection defined in Fig. [T] In such a case one could find 
novel envelope functions ^k,m{y, z) associated with every 
quantum state, with the following symmetry for E\/2 or 
2 E 1 / 2 states respectively : 



(-i) i+m <m /2 (y>-*) 

(-l)i+-+l^(y,- 2 ) (g) 



To illustrate this we display in Fig. [4] the contour plots for 
the same ground state as in Fig. [3] Although seemingly 
different, this new envelope function representation for 
the eigenstate carries exactly the same physics, i.e. gives 
the same expectation values for all physical observables. 
It should be mentioned that the basic reason for the m- 
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Figure 4: Contour plots of valence band envelope functions 

^ k ^ ( y, z) of the ground state of symmetry 1 E 1 /2 in a typical 
V-shaped QWR, when the Bloch function basis diagonalizes 
the J z component corresponding to the [110] crystal direction 
(same QWR as in Fig. H]). (a-d) m = 3/2,1/2,-1/2,-3/2 
function components respectively (J z ). 

behaviour of the envelope functions in Eqs. © compared 
to Eqs. l[7]). which might seem surprising at first, can in 
fact be explained in a very intuitive way by looking at the 
behaviour of angular momentum components through a 
planar reflection a : indeed the sign of the in-plane com- 
ponents of the angular momentum are reversed, while the 
perpendicular component is conserved, i.e. 



erJ(7 



-J + 2 J z e 2 



(9) 



Therefore it is obvious that if one uses a Bloch basis 
diagonalizing the component of the pseudo-angular mo- 
mentum J perpendicular to the symmetry plane (i.e. J z ), 
every envelope function will be mapped onto itself (either 
in a symmetric way or antisymmetric way), whilst the 
+m and — m components will be mapped onto each other 
if one diagonalizes J y . To show that the envelope func- 
tion spinors linked with the two double group irreps 1 E 1 / 2 
and 2 Ei / 2 have opposite alternating parity in Eqs. © re- 
quires a more detailed analysis. However we prefer to re- 
legate this discussion after the presentation of the general 
theory, since most properties will become obvious. 



2. From low to high symmetries : C<xv, Cs v and higher 

We have thus shown that a careful choice of basis al- 
lows, in the case of C s , to symmetrize individual envelope 
functions. Would this be possible in the case of higher 
symmetry? We intend now to clearly demonstrate that 
the approach suggested in Ref . [l(| , where we introduced 
the concept of Optimal Quantization Axis (OQA), has 
limits, motivating the more elaborate MSRF approach. 
We shall now take quantum dots as simple examples. 
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Let us start with C 2v , the next higher symmetry de- 
picted on Fig. [SJ A practical example is pyramidal InAs 
QDs which are have the shape of a rhombus-based py- 
ramid. We shall denote the two perpendicular symmetry 
planes o~ y and o~ z for compatibility with the axes used 
in this paper. Since C s is a sub-group of C211, we could 
use the same basis diagonalizing <r z , and reduce the pro- 
blem to the half-domain. However it would be highly de- 
sirable to use the additional symmetry with respect to 
a y to further split the domain of solution. No problems 





Figure 5: Perspective view of a G^v QD 

for electrons, but problems would arise for holes, since 
a 4 band k ■ p model would be required, and from our 
previous discussion of C s symmetry we immediately see 
that if we would take the spin quantization axis along 
z, the spinors envelope functions would become alterna- 
tingly even/odd with respect to z, but with respect to 
y would necessarily obey a symmetry relation coupling 
+m and —m (c.f. Eq. (JTJ). Oppositely if one would have 
chosen the y basis, the symmetry relations with respect 
to z would have become badly behaved. Therefore it is 
apparently never possible to obtain symmetric envelope 
functions in the two directions simultaneously, and solve 
on the half domain in the two directions ! 

There is a rather simple explanation at a more fun- 
damental level for this different behaviour. For the elec- 
tron symmetry, described by the single group, a y and a z 
commute, and therefore one can in principle diagonalize 
the two operations simultaneously, and get simultaneous 
good quantum numbers linked with them. For the hole 
spinorial symmetry, described by the double group, a y 
and a z do not commute, indeed the general commutator 
can be written as : 



[o-y,a z }=C 2 [l + {-l) 



(10) 



showing that when j is half-integer (here j — 3/2), it is 
never possible to diagonalize simultaneously both sym- 
metry operations (such a fact can also be related to the 
appearance of a 2D irrep for the double group and the 
properties of its corresponding 2x2 unitary matrix repre- 
sentation D El / 2 (fi),g £ C^v). As a result of this analysis 
we suggested in [10( that the optimal basis was naturally 
the one diagonalizing the projection of angular momen- 
tum along the third perpendicular axis, allowing to treat 
a y and o~ z on an equal footing, and diagonalizing the ro- 
tation C 2 — G y ■ o~ z . However a closer inspection reveals 
that no solution on the quarter of the domain is yet al- 
lowed in C 2v by this idea, this will only become possible 
with MSRF. 



Even more challenging is the next higher symmetry : 
Cs v symmetry, with three symmetry planes like in Fig. [5] 
Such a symmetry is also of practical interest [13, EH • The 




Figure 6: Perspective view of a C$ v QD 

reference axes of the crystal, and our x, y, z labels, are 
shown in more detail in Fig. together with the three 
vertical symmetry planes o~ V i,i = 1,2,3. In addition to 

i y [112] 
1 1 



Ov3 



[110] 1 <- 



Gv2 



© X [111] 



* * A* 

I 
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Figure 7: Axes and cross section of the Al x Ga\~ x As Cz v QD 

the three improper rotations (vertical symmetry planes) , 
called o~ V i , i = 1,2,3, the C^ v group includes as addi- 
tional symmetry operations two rotation of ±120° (C3 ) 
and the identity (E). 

The Cq v group displays two ID irreps A\ and A 2 , and 
one 2D irrep (E), even for the single group, leading to two 
basis function (partner functions) for the subspace rela- 
ted to degenerate eigenvalues. There is a supplementary 
difficulty linked with this 2D irrep, in particular the cor- 
responding 2D matrix representation explicitly depend 
on the basis functions. The simpler ID single group ir- 
reps A\ and A 2 are respectively even and odd with res- 
pect to all the symmetry planes. Therefore, for electrons 
it is straightforward to see that one can compute easily 
A\ and A 2 eigenstates via the solutions on 1/6 of the 
domain by imposing respectively Neumann or Dirichlct 
boundary conditions on the two symmetry planes. For 
the degenerate E irrep one can choose the two basis func- 
tion such that they are either even or odd with respect to 
one of the symmetry plane (mirror), for instance a~ v \, but 
it will not be possible to diagonalize simultaneously any 
two of the mirrors at the same time, for instance a v ± and 
a v2 , since they do not commute, as can be seen in the 
multiplication table (see annexe fA")) . Therefore it is not 
obvious how to use the mirror symmetries to solve on a 
reduced domain smaller than one half. One may conclude 
that in a C% v heterostructure the problem that appeared 
only for valence-band holes in the case of C 2v symmetry 
already appears for spinless electrons : it is not possible 
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to solve with this technique on the most reduced domain, 
which, in C$ v , should be smaller than one half of the full 
domain ! For holes in C 3v there is also a 2D faithful self- 
conjugated double group irrep (£1/2)) but also two ID 
mutually conjugated irreps { l E 3 / 2 , i = 1,2) [2l| . 

Finally all these questions become more severe in hi- 
gher symmetry, like the hexagonal Cq v symmetry (c.f. 
perspective view of a C§ v QD in Fig. [5J which is also of 
practical interest [l9|, H3] ■ A characteristic of this group is 
that the double group displays only 2D degenerate irreps. 
Although such structures have been discussed to some 
extent in the literature with the help of group theory [T|| , 
the symmetry properties of the envelope functions have 
never been studied and discussed. With MSR it is pos- 
sible to approach systematically all higher symmetries, 
and for Cq v the results will be given in subsection I VIII Al 




Figure 8: Perspective view of a Cev QD 

To summarize, we can identify a limit between low and 
high symmetry groups in our sense : the appearance of 
a 2D irrep, which is a manifestation of the non-abelian 
character of groups like the dihedral groups (C nv ) that 
we have considered. Such irreps with dimension greater 
than 1 complicate a lot the question of the symmetry of 
the basis functions, making it non-trivial. 

Finally, the concept of spin quantization axis direction, 
linked with optimal pure rotations of the original spino- 
rial basis, is not a suitable concept to tackle such higher 
symmetries (with spin, already Civ 0- Clearly maximal 
symmetrization of the envelope functions could not be 
achieved, and computation on a reduced domain was not 
enabled. In the following, a radically new approach for 
HSH is presented which will fulfill these goals. The opti- 
mal spinorial basis is obtained by a more general unitary 
transformation corresponding to the reduction (block- 
diagonalization) of the spinorial representation and rela- 
ted to the choice of double group labelled basis functions. 

In the next chapter we shall start developing the MSRF 
formalism from beginning, by looking at transformation 
laws in "orbital" and "spin" space. Whenever needed the 
general theory will be illustrated by the typical case of a 
Cs v QWR (2D problem), either for the spinless conduc- 
tion band, or for the valence band, again with the four- 
band Luttinger Hamiltonian. Numerical examples were 
worked out in real space with a FE approach incorpora- 
ting linear elements only. 



III. GENERAL TRANSFORMATION LAWS, IN 
ORDINARY AND SPINORIAL SPACE 

The study of transformation laws under symmetry ope- 
rations is a prerequisite for the efficient use of symmetry 
and group theory explicitely on a given problem. One 
needs to know exactly how conduction and valence band 
envelope functions do transform under symmetry opera- 
tions. This is of course related to the corresponding k ■ p 
Hamiltonians presented in Sec. [TXJ In this section we in- 
troduce in details transformation laws, which are also a 
prerequisite for the development of our new theory, in- 
dependently of any symmetry consideration. For clarity 
we treat seperately scalar functions and spinors. In the 
last part of the section only we introduce symmetry, and 
look at the resulting constraints on envelope functions 
encountered in k ■ p theory. 

A. Transformation laws 

Let us first introduce transformation laws for simple 
spinless scalar functions (typically the quantum wave- 
function of an electron in a conduction band) , and then in 
a subsequent step expose the transformation laws in the 
spinorial case (typically a hole in the valence band). To 
eliminate any ambiguity in the following we shall always 
use a passive point of view for the symmetry operations, 
i.e. the operations are always considered as coordinate 
transformations linked with a change in reference frame, 
i.e. they are not rotations of the physical system. 

1. Coordinate transformations 

Let us first define a basic transformation of coordinates 
c linked with a change in orthonormal cartesian reference 
frame. It is defined by an orthogonal 3x3 matrix 5ft (c) 
belonging to 0(3) and defining the basis vectors of the 
new frame with respect to the old one : 

where i,j € {x, y, z}. It is either a rotation or an im- 
proper rotation, any rotation can be parametrized by its 
Euler angles a, /? and 7, and the corresponding set of 
rotation matrices 5ft (g) = 5ft (a (3 7) define a representa- 
tion of the rotation group (S*0(3)). These matrices can 
be systematically constructed using the generators of the 
rotations [32[ : 

5ft (a /3 7) = e +t ~< J *e +if}J ye +2aJ * (12) 

where J Xl J y , J z are the components of the angular mo- 
mentum pseudo- vector J. Improper rotations can always 
be decomposed as the product of the spatial inversion i 
with a proper rotation, therefore their matrix represen- 
tation is written as the product of 5ft (i) — —I3 with the 
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corresponding 5ft (a f3 7) as defined in Eq. (|12l A typical 
example is ag, a mirror symmetry with normal s, which 
is the product of the inversion i with C2(s), a 7r-rotation 
around the axis s . Therefore 5ft (<7j) — — 5f (7T, s), where 
temporarily the Euler angle notation is left out. Simi- 
larly, in the following we shall use c to denote the change 
in coordinates linked with any arbitrary change of ortho- 
normal reference frame, and often 5ft = 5ft(c) (without ar- 
gument) will denote implicitely the corresponding ortho- 
gonal matrix representation (the argument of 5ft(c) will 
be restored only when absolutely needed for clarity). 

Let us now recall how the set of components r = 
(x,y,z) of a vector r with respect to a given basis 
{e{\ ,i G {x,y,z} transform under a change in coordi- 
nates (new basis {e^}) : 



guments (r, k) should thus be understood as follows : 



XtijXj ^ r' ftr 



(13) 



which is a contragredient law with respect to ea. (|ll|) . as 
it should for a passive point of view. 



2. Transformation of scalar functions 

In the Hilbert space rl — C 2 (R d ) corresponding to the 
set of possible electronic wavefunctions in d confined di- 
mensions (d=l,2,3), one can associate linear operators d c 
to every possible coordinate transformation c : 



d c : rl — > rl 



(14) 



The new mathematical function ip' is a function of the 
new coordinates, but is defined as representing the same 
quantum state, i.e. 

V/(r') = ^(r) => ^V) = $c m (!•') = V(« _ V) (15) 

It is easy to check that this definition leads to the expec- 
ted multiplication rule i? C2 ci W\ ( r ) = $c 2 &ci [ip] ( r ) = 
§ C2 ^o5ft-!( Cl )] (r) = ^- 1 {c l )^.- 1 {c 2 ))v) = 
V>((5ft _1 (c2Ci))r) for two successive coordinate transfor- 
mations Ci and C2. 

Let us now consider the generic form i/(r,k) of the 
scalar k ■ p Hamiltonian in mixed position and momen- 
tum representation (of complementary dimensions d and 
df = 3 — d respectively) which appeared in Section III Al 
(Eq. {I)). This generic form is applicable for heterostruc- 
tures of any dimensionality df, from 3D to OD (df is the 
dimensionality in fc-space related to the "free-like" mo- 
tion of charge carriers in the directions with full transla- 
tionnal invariance at the heterostructure level). The ar- 



bulk (3D, d=0) :(r,k) 
quantum well (2D, d=l) :(r, k) 
quantum wire (ID, d=2) :(r, k) 
quantum dot (OD, d=3) :(r, k) 



(k — [k x , ky, k z )) 
(z,ky = (k x ,k y )) 
(r ± = (x,y),k\\) 
(r = (x,y,z)) 

(16) 



The new transformed Hamiltonian operator is obtained 
by enforcing (in Dirac notation) H \ip) = H' \ijj'), 
which leads to 

H'(r', k') = d c H(r', k') [i? c ] _1 = (5fT V, 5fT x k') 

(17) 

Note that k is implicitely understood as the vectorial 
components of the wavevector (covector). In Eq. (jTTJ) the 
operations c are always considered in 3D (i.e. belonging to 
0(3)), therefore the eigenfunctions of the generic Hamil- 
tonian must have an auxiliary index in order to transform 
consistently, i.e. 

^ (r') = d c [W] (J"') = ifa-xv (3T V) (18) 

Such a generic scalar Hamiltonian is typically a quadra- 
tic form of the momentum, i.e. of the components of 
ky along the non-confined directions (translational in- 
variance), and of differential operators which are kept 
in the confined directions (according to the correspon- 
dance kj^ = — iVj_). The parameters of this Hamilto- 
nian, i.e. the effective masses and the confinement po- 
tential, are functions of positions like in Eq. Jl} . It may 
also include anisotropic masses as follows : i/(r,k) = 
Hq(y) + k 4 C(r)k, where C(r) is the "matrix" (tensor) 
of coefficients and Ho(r) is a scalar operator. In this ra- 
ther general case the new Hamiltonian is obtained by 

H'(r, k) = H^-h, 5JT *k) = Hq^t) + k* C'(r) k 

(19) 



where 



C'(r) = 5RC(5ft~ 1 r) 5JT 



(20) 



Is easy to understand equation (12"0"1) when C is inde- 
pendent of r : k* C k is a scalar, invariant under a pas- 
sive transformation and k are the vectorial components 
of wave vector, then C is a tensor two times covariant 
and, for this kind of tensors, the tranformation laws give 



Ers^rC rS 5ft 



3. Transformation of spinors and spinorial fields 

In a model including the spin of the charge carriers, 
the quantum state may generally be described by a spi- 
norial field ip(r). In this section, for clarity, we will use 
the underscore to explicitely denote the spinorial charac- 
ter, which is assumed to be of dimension 2j + 1 where j 
is a strictly positive half-integer. In the k ■ p multiband 
envelope function formalism, the spinorial components 
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are related to the Bloch function basis at a high sym- 
metry point of the Brillouin zone, diagonalizing the Ha- 
miltonian of the charge carrier in the bulk crystal struc- 
ture [33l ] . 

The corresponding Hilbert space is spanned by the tensor 
products of spinors belonging to C^ 2j+1 ^ with the enve- 
lope functions belonging to £ 2 (R d ). Hence when consi- 
dering truly arbitrary coordinate transformations one 
should consider to perform the transformations indepen- 
dently in both spaces. 

Let us first look solely to SU(2) spinor transforma- 
tions. By contrast to vectors in normal cartesian 3D 
space (c.f. Eq. (flTt ). spinors obey different transforma- 
tion rules under rotations, called spinorial transformation 
rules [38|, HH . For pure rotations the transformation rules 
can always be related by a similarity transformation to a 
set of matrices (a /3 7), indexed by the Euler angles, 
and called the Wigner representation |32l EbT ] , which is a 
(2j + l)-dimensional projective representation of SO (3) : 



W U) (a/3 7 ) = e +n ± 



«) +iaJ (. 
e —ve =» 



(21) 



When j is half-integer it is a spinor representation of 
SU(2) (instead of 5*0(3)). The double underscore no- 
tation introduced here always denotes a square matrix 
character in spinorial space. A typical feature of spi- 
nor representations is that a 2tt rotation around any 
axis s will always be associated to a sign change, i.e. 

(2ir, s) = — 1. For the representation of improper 
rotations one can again use the factorization of the in- 
version, however a small complication is the Wigner re- 
presentation of the inversion, which can be either ±1. We 
shall not discuss here the proper choice of the latter sign 
since it will not matter within this paper. In some cases 
it might be important, see |28j |. 

A crucial feature of our approach is now to consider 
more general unitary change in coordinates c, not linked 
with 3D rotations, characterized by the general set of 
square matrices V_(c) of dimension 2j + 1 belonging to 
the group U(2j + 1). 

Let us start with the (2j + l)-dimensional Bloch functions 
basis denoted {\j, m)}, as is customary in the field, and 
a new basis m)'} differing by a rotation c = (a/37), 
so that : 



\j, m ) 



(22) 



which is similar to the standard transformation law of 
partner functions linked with the irrep of dimension 2j+l 
of the rotation group under c _1 . It is important to make 
here two remarks for clarity. First we stress that for such 
Bloch functions j simply correspond to a label related to 
the transformation law (f2"2"j) and to the dimensionality of 
the Bloch functions basis, but in this context it is not a 
true angular momentum quantum number (it just refers 
to the transformation properties, and not truly to rota- 
tions). A second important remark is that the {\j, m)} ba- 
sis, and its corresponding transformation law under rota- 
tions (|22p . already incorporate, in addition to pure spin, 



an orbital part, since the "bulk" spin-orbit interaction is 
already diagonalized by the Bloch function basis. Later 
below a purely orbital transformation, operating solely 
on the envelope function, will have to be simultaneously 
added. When one considers more general transformations 
of coordinates in spinorial space (c £ U(2j + 1)) : 



A3) 



ip 1 — > ip 



. c(2i+i) 



(23) 



the spinor components ip' can be transformed using the 
matrix equation 



tp' = V(c)tp 



(24) 



which is anologous to Eq. (fT3"|) . The matrix V_(c) defines 
a more general basis change than the one specifed by 
Eq. (EU). 

Let us now consider all the possible transformations of 
spinorial fields, which are defined by two separate coor- 
dinates transformations in both spaces c\ £ 0(3) in real 
space and C2 £ U(2j + 1) and combine them in d c as- 
sociated with c = (c\,C2) acting in the tensor product 
Hilbert space : 



(25) 



where #c? D ' ) is the operator defined by Eq. (fT4")) and $ 



U) 

C2 

by Eq. (|23| . Here we have kept the possibility of arbi- 
trary and different changes of coordinates in both spaces, 
which will prove crucial later. Under a coordinate trans- 
formation i? c an arbitrary quantum state of our system, 
described by a spinorial field, transforms therefore like 



21k v 



H C 2)^- 1(ci)k 



(» _1 (ci)r) 



(26) 



We see clearly in Eq. (|26|) that the spinorial character of 
the transformation docs couple different envelope func- 
tion components through the V_ matrix, and is more com- 
plicated than the simple 3D transformation of the enve- 
lope functions appearing in Eq. (fTSJ). This complication 
is one of the roots of the difficulty in interpreting the 
individual symmetry of envelope functions in a heteros- 
truture with a given symmetry. 

In the spinorial case one can transform the Hamiltonian 
operator matrix too, but one must take into account the 
presence of the V_ matrix in a similar way : 

£'(r,k) = tf c il(r,k) [tfep 1 

= ZEj&r^tirtyYr 1 (27) 

where, in the last two lines, we have used again the shor- 
thand notation = 3?(ci) and V_ = V_(c2)- 
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B. Symmetry and resulting constraints on 
eigenstate envelope functions 

Let us now assume a heterostructure with a given 
symmetry group Q = {g} of spatial transformations. 
More precisely, if there is translationnal invariance (i.e. 
df > 0), let us take Q as the small point group of k = fey 
of the structure, with cardinality \Q\, and defined by the 
restriction : 



k = k 



(28) 



The operations g € Q are either pure rotations or roto- 
inversions (mirrors). The heart of the passive point of 
view is just to express than any full coordinate trans- 
formation corresponding to a symmetry element g of the 
structure will leave invariant the form of the k-restriction 
of the conduction and valence band Hamiltonians (hence 
from now on we shall leave out the implicit k-subscript on 
eigenstates and envelope functions appearing in Eqs. (|18[) 
and (I26|) . 

For the conduction band the product g" of two symmetry 
operations g and g' will simply follow the multiplication 
table of the single point group G, whilst for the valence 
band one must necessarily use a double group notation 
of symmetry operations g £ Q due to its spinorial nature 
(please note that here we use the notation g to denote any 
element of Q 1 instead of the single group element g multi- 
plied by a 27r-rotation, as is quite standard (see e.g. [2l|)). 
The composite index c for the corresponding full coordi- 
nate transformation can be identified with g, provided 
one understands it as {g,g) where g is the single point 
group image of g. Clearly this composite index obeys also 
the double point group multiplication table since the pro- 
duct is defined as (g,g) ■ (g',g') = {g-g',g.g') = (g",g")- 
Using Eq. (|27| let us now express the invariance of the k- 
restricted Hamiltonian with respect to a given symmetry 
operation g of the coordinate system : 

M!(r, k) = VJf (8T x r, k) VT 1 = K(r, k) (29) 



where shorthand notations have been used. Eq. (|29|) is 
equivalent to state that every symmetry operation g com- 
mutes with the k-restricted Hamiltonian, which allows 
the use of a well known theorem [33| that states that 
every eigenspace of the Hamiltonian can then be labeled 
by an irreducible representation (irrep) r of Q (meaning 
that except for accidental degeneracies its dimension df 
is necessarily the dimension of the irrep) , and that a basis 

of partner eigenstates ip T (r) (/x = 1, . . . , df is called the 
partner function index) can be found such that under a 
basis change dg it transforms according to 



d 



(30) 



where the set |_D r (g)| form a unitary irreducible matrix 
representation of the irrep T. 



From the physical point of view, Eq. (|30| means that a 
transformed symmetrized eigenstate V^( r ) under $g does 
not only have the same energy but can also be developed 
on its partners, so that there are constraints on indivi- 
dually transformed envelope functions since they must 
"reconnect" on all the other ones in a very intricate fa- 
shion. Using Eq. (|2l)|) we find 



]TLd F (.9)1 £(t) = V(9)€(^ 1 (~9>) ( 31 ) 

L — ' L J v^i — " — — f 1 



This equation forms the starting point of our theory. 
It will be exploited in the next section where we use 
the explicit separation of the spatial and the spinorial 
part of the operators in Eq. (|25p . and then seeks the 
basis that will minimally reconnect envelope functions 
according to (|3"Tj) . 



IV. THE FULLY SYMMETRIZED OBB BASIS 
AND THE SEPARATION OF SPINORIAL AND 
SPATIAL PARTS 



By definition, a 3D representation of any point group 
can be obtained simply from the analytical expression of 
the 5i(a, f3, 7) matrices (|T2"j) and by factoring out the in- 
version for improper operations. In spinorial space things 
are slightly more complicated since, depending on the ba- 
sis used to express the Hamiltonian, the V_(g) matrices 
do form a spinorial (2j + l)-dimensional representation of 
the double group Q, but not necessarily given by Wigner 
matrices (|2T1). 



In this section we shall start from a spinorial Bloch 
function basis which, at high symmetry points of the 
Brillouin zone, is often denoted m)}, due to its 
symmmetry transformation properties. For the valence- 
band of semiconductors like GaAs, at the so-called 
T-point, one restricts to j = 3/2 [22j giving rise to the 
Luttinger Hamiltonian (c.f. Eqs. © and ((H)). However 
the formalism holds for an arbitrary large number of 
bands. The heterostructure symmetry is kept general, 
but illustrated with C^ v symmetry. A few other cases, in 
particular the Cq v , D^h , C n , C s symmetry groups, will 
be shortly discussed in Sec. IVIII1 

To introduce the concept of the Optimal Bloch function 
Basis (OBB) with respect to a heterostructure with a 
given symmetry, we must heavily rely on the explicit 
separation of orbital and spinorial part carried out in the 
last section, and on the possibility of separate coordinate 
transformations in both spaces. 
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A. The optimal Bloch function basis 



Our main goal is to simplify Eq. (1311) and minimize 
the coupling between different envelope functions. The 
main novel idea is to perform once for all a unitary co- 
ordinate transformation c = (l,c 2 ),c 2 G U(2j + 1), i.e. 
purely in spinorial space, corresponding to a best choice 
of the Bloch function basis, such that the set of reducible 
matrices V_(g) appearing in (|3"Tj) would become block- 
diagonal ! 

Let us first note that some aspects of this idea are 
not completely new. Up to now one as used 3D rota- 
tions of the quantization axis (direction of J z ), para- 
metrized by the Euler angles c 2 = (a, ft, 7) such that 
J z would be transformed towards J' z , diagonalized by 
\j, m)' [13, [HI- The corresponding image of e z is the so- 
called Optimal Quantization Axis (OQA) [l(J. The more 
trivial case of quantum wells grown in [hhk] direction, 
where the OQA is always [hhk], is treated in |31|. The 
OQA was supposedly considered the best choice [T(| to 
simplify the Luttinger Hamiltonian. Indeed a change in 
basis \ j, m)' — ^2 n U nm (c^ ) \ j, n) has lead to a new form 
of the Hamiltonian expressed in the new basis as 



H'(r,k)=g( C2 )ff(r,k) [g(c 2 )] 1 



(32) 



which lead to simpler wavefunctions in the C s case. It 
also lead to a corresponding new representation of the 
symmetry operations 

1 -1 



K'tf) =££(<*) Z(s) [£(<*)]" 



(33) 



The optimal choice of angles c 2 = (a (3 7) was carefully 
made [13, HH, and discussed partly in Sec. [TTJ but 
we have seen that it cannot be efficient for higher 
symmetries. 

In the present approach, the novelty is to consider more 
general unitary transformations c that cannot be repre- 
sented as rotations of the original reference frame. For 
this purpose we shall seek the coordinate transformation 
c = (l,c 2 ),c 2 G U(2j + 1), but possibly c 2 <£ SU(2), to- 



wards new Bloch states 



T,fiJ, where T is the label of a 

double group irrep, and /i a corresponding partner func- 
tion label. This set of states, that we shall call Optimal 
Bloch function Basis (OBB), is a fully symmetrized Bloch 
basis, i.e. symmetrized according to the symmetry of the 
quantum hetero •structure. From its definition 



r6^) = ^f/ m;F ^(c 2 - 1 )|j, 



(34) 



The notation makes it clear that every new basis state 

Tb,p\ will map under passive symmetry operations like 

the standard set of partner functions of an irrep of the 
double group : 

dp 



n 



U) 



S'=l 



/3'0 



(35) 



but with g~ l on the right hand side. Let us now look at 
the transformation properties of the new spinorial com- 
ponents of the field tjj B (r) = i?c'0( r ) — IZfe) ^(r), trans- 
forming as 



^V B (r) 



V B (g)r^- l (9)r) 



(36) 



From Eq. (|35|) it is clear that the change in basis speci- 
fied by the t/(c 2 ) matrix induces automatically a new 
block-diagonal representation of the symmetry opera- 
tions V_ B (g) where Eq. (|3"3"|) then becomes 



(37) 



The ns = J2b n f b blocks of dimension dp of V_ (<?) are 
labelled by Tf, and a possible multiplicity index runnig 
from 1 to n^ b . Indeed a given representation Tf, may ap- 
pear more than once (a pratical example is the C s group 
rediscussed later), however for simplicity we shall forget 
it since it can be restored without difficulty. 
Let us now separate every irreducible blocks Tf, of V_ B (<?), 

and of every partner eigenstate ip B (r) of the irrep T : 

— A* 



if (5) = Ef ® 



(38) 
(39) 



It is understood here that every V^ b (g) has only a single 
non-zero square block on its diagonal, which we shall de- 
note V Tb (g) = D Th (g) to pinpoint that its dimension is 

d^ b instead of (2j + 1). Similarly ifj r ^ rb (r) is defined as 
having only a single subset of tip relevant components 
which we shall denote ip^ l ' Th (r) (all other components 

are zero). If one now applies the transformation c 2 to 
Eq. ([3Tj) . and uses Eq. ([37| . and identifies every subblock 
labelled by Tb, one finds : 



E 

i/=i 



tu' Li (r) = Y L "(g)^ L "(U-\g)r) (40) 



Clearly this equation shows that our goal to minimize 
absolutely the coupling between different spinorial com- 
ponents under symmetry operations is achieved : the cou- 
pled components are reduced and grouped according to 
the irreps r&. However this is not sufficient, and more 
work needs to be done, and will be carried out in subsec- 
tion |lVCl to ensure that the resulting envelope functions 
have maximum symmetry. Let us first illustrate the de- 
velopment of this section with a concrete example. 
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B. The OBB for the Luttinger Hamiltonian in Cz v 
symmetry 

The matrix representation of point group symmetry 
operations for the standard 4x4 Luttinger problem at 
hand are in fact naturally the original Wigner repre- 
sentation 32 1, i.e. V(g) = Wy> (a ftj), since they are 



consistent with the J matrices appearing in ([2]). One 
just needs the reduction in block diagonal form for C 3v 
symmetry which can be, in this case, efficiently carried 
out with the help of the matrix traces. One then finds 
that this matrix representation of C 3v is reducible to 
1 £ ; 3/2 © E1/2 © 2 E 3 / 2 - The ad- hoc Z7(c 2 ) transformation 
can be found with standard techniques [21| and leads to 
the following fully symmetrized Bloch function basis 



E 



1/2 



E 



1/2 



3/2) 
2) 



^3/2) 



¥ 

1 

2i' 
1 

— t 

2i 



3 3 
2' 2 
3 1 
2' 2 
3 

2'" 
3 

2'' 



V3 



3 

2' 

3 

2' 



V3 



1 

"2 

3 
2 
3 3 
2' 2 
3 1 
2' 2 



,(41) 



where the \ j, m) basis pertains to the a;, y, z axes in Fig. [7J 
The reduced block-diagonal form of the V_ B (g) matrices 
is then : 



V B (9) 



x E3/2 {g) 











x E3/2 (~g) 



(42) 



Here X E3/2 (sOj i = 1, 2 are the characters of the ID irrep 
^3/2) i = 1)2 respectively, and D El / 2 (g) is a 2D matrix 
representation for -E1/2 that one can freely fix by choosing 
the two partner functions. We have specifically chosen 
also 



D^(g)=D*(g) X ^(~g) 



(43) 



where 
group 
sentations 



D E (g) is the 2D representation for the single 
irrep E, according to the product of repre- 



E ® 2 E 3 / 2 ~ Eif 2 (see appendix [A)- It is 
important to realize at this point that there is a degree 
of arbitrariness that cannot be avoided : for example the 
order of appearance of the blocks in (|42"T) , or the fact that 
both matrix representations D E (g) and D El / 2 (g) can in 
principle be chosen independently. To justify our specific 
choice two remarks are in order : first it is desirable 
to keep the p, q, r, s-form of the Luttinger Hamiltonian 
(linked with the form of the time- reversal operator), and 
this can be achieved only with a few possible orderings 

of basis states Tb,ft^, second Eq. (|43|) is a particular 

restriction whose motivation will appear more clearly 
later. 

An alternative group-theoretical view on the decompo- 
sition of the Wigner representation into block-diagonal 



form (|42|). is to consider the diamond group Oh, which 
is the symmetry group of the Luttinger Hamiltonian 
linked with the underlying bulk semiconductor. With 
respect to the symmetry Oh, the C 3v heterostructure 
can be viewed as a symmetry breaking perturbation (due 
to the mesoscopic heterostructure potential), therefore a 
good basis within the 4D subspace linked with the irrep 
Es/2,g °f Oh will be found using subduction tables [211 ] 
which give F 3 / 2 . g -> ^3/2 © E 1/2 © 2 E 3/2 . 

It is interesting to note that the preliminary co- 
ordinate transformation c = (1,62) is here purely 
in spinorial space, but in fact this is not a general 
feature. One can explain it in the following way : for 
C 3v the 3D matrix representation {3?(g)} is already 
in a reduced block-diagonal form (A\ © E) with the 
basis presented in Fig. [7] (e x is invariant respect to 
every symmetry operations of the group and e y ,e z 
are mutually coupled according to E irrep). In fact 
here one has directly implicitly chosen the vectorial 
basis according to the symmetry of the heterostructure 
analogously to the Bloch function basis in Eq. (|42|) ! A 
counterexample would be in the C n group, which will 
be shortly discussed in section I VIII CI : one indeed needs 
to introduce a slightly more complex preliminary coor- 
dinate transformation of the form c = (01,02) such that 
the matrix representations {3?(g)} and (H(g)} would 
become simultaneously block-diagonalized according to 
the irreps of C n . 



Ultimately reduced envelope functions in the 
OBB basis 



We now come to the last positive by-product of the 
introduction of the unitary coordinate transformation c 
towards the OBB : the ability to define "ultimately" re- 
duced envelope functions (UREF). 

Let us revert Eq. (|4Tfl) by multiplying it by V^ b (g~ 1 ). 
We now look at each individual envelope function com- 

i.r.f, 



ponents V' p '^ i> ( r ) of ip^ ' 6 (r) : 



,r,r ( 



(^ 1 (3)r) = 



EE 

0>=1 u=l 



[D f (g) 



(44) 

where the envelope function components are now clearly 
related to the OBB basis (j34|) through the indices Tb and 
ft. One should now make a fundamental remark, namely 
that Eq. (|44j) can be interpreted as a reducible transfor- 
mation law for individual envelope function components 
under the symmetry group of the heterostructure ! In- 
deed the transformation matrix clement can be written 



DIM D^^g) 



which contains all elements of the 
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product of representations 

r®fj = ; 



T,r; ; r, 



(45) 



where n 



r,r; ; r 



is the multiplicity of the T a irrep in the 
product. Moreover, since T a appears in a tensor product 
of the double group representations, r a is necessarily a 
single group representation ! (note also that for point 
groups p*. r < 2, however we shall deliberately ignore 
in the following the corresponding additionnal multipli- 
city index : first it is trivial to restore it if needed, se- 
cond one does not usually need it (simple-reducible point 
groups). 

Let us therefore introduce unitary generalized Clebsch- 

Gordan coefficients C^'^^" which perform in practice 
the block decomposition of the reducible matrix appea- 
ring in (14"4"1) . i.e. 



E 



E 

i fi' i 



c 



(46) 



This equation naturally leads us to introduce, by the use 
of Eg. (1451) . a set of Ultimately Reduced Envelope Function 



components (UREFs), denoted c/>~' ra (r), and associated 
with every subspace of the nanostructure Hamiltonian of 
symmetry T and every block Tb ■ 



&)=EE«r° 



(47) 



M=l/3=1 



This equation can be reversed easily by using the uni- 
tarity of the Clebsch-Gordan matrix, which gives a de- 
velopment of the original envelope functions in terms of 
UREF's : 



E [ C M,/3;« r " 



4: r :( r ) 



r b ,a 



(48) 



Fourth, one interest of Eq. (|48| is to suggest an alter- 
native way to arrive at UREF's by imposing for every 
eigenstate a fixed variance development on the fully sym- 
metrized OBB (a generalization of the invariant deve- 
lopment approach for Hamiltonians [ll|). Fifth and last, 
the definition of the Clebsch-Gordan matrix entails ad- 
ditionnal phase factors due to the conjugation of Tf, (c.f. 
Eq. gSD). 



D. The UREF's in C 3 „ symmetry 

We now illustrate the general results of the preceding 
section by studying again the specific case of a C^v hete- 
rostructure. Let us consider hrst the non-degenerate irrep 
r = 2 E 3 / 2 . The first component of ip Ea/2 (r), related to 
the 1^3/2) basis function, reduces directly, via Eq. (|44l) . 
to the simple uncoupled expression : 



X A2 (2)V<r /2 (r) 



(50) 



where we have used the fact that z E l 3/2, i = 1, 2 are mu- 
tually conjugated irreps, and that the direct product re- 
presentation 2 i?3/2 ® 1 ®3/2 = ^2 (c.f. table in appen- 
dix [X]) . Therefore we finally obtain, in agreement with 

2 E 

Eq. (|4"9")l that ip-. 3/2 (r) transforms like the single group 
irrep Ai. In a completely similar way, it is possible to 

2 E 

show that the last component of ip 3/2 (r) transforms like 
Ai and that the two central functions like the 2D mutual 
partner functions of the irrep E. Therefore one can write 

2 E 

the eigenstate ip 3/2 (r) in most simple symmetrized form 
as : 



ip E3 / 2 (r) 



( 



<^(r)\ 
Al (r)/ 



(51) 



The definition of UREF's, Eqs. fl47|) and Q48J, deserve a 
number of important comments. First, the UREF decom- 
position (|47p is totally general. Second, for every subspace 
of the Hamiltonian related to the symmetry T there are 
rip = J^r f 6 n r b n f f*T se ^ s °^ UREFs, in particular rip 
can be higher than (2j + 1), note also that the UREFs 
are characterizing strictly T and r?,, i.e. they are inde- 
pendent of the partner function indices /x and /?. Third, 
from Eq. (|45l) one can show that every set of UREFs 
indeed transform minimally under symmetry operations, 
simply like the partner functions of the single group ir- 
reps T a : 



r b ,a 



(r) = ^m^Cg)r) 



r b ,a 



(49) 



where we have left out for clarity the Ft, label on the cp 
functions (it can be easily restored by using the conse- 
cutive labels of the Bloch function basis). For the £3/2 
irrep we obtain a similar decomposition : 



tp B 3/ 2 (r) 



V 



\ 



0f(r) 
-<t> A2 {v) ) 



(52) 



It should be pointed out, that the <p- functions appea- 
ring in Eq. (|52[) are not identical to those appearing in 
Eq. (f5"Tj) . this becomes obvious if one restores the implicit 
r label. However time-reversal will induce a mapping bet- 
ween these functions (see e.g. [Hj]), but simultaneously 
with a change of the k- index to — k (we have not used 
here the complicated notations of the preceding section 



15 



for the (^-functions because they are unnecessary in this 
fairly trivial case) . For this reason we also included a mi- 
nus sign in (f> A2 (r) in Eq. (f52j) . 

For the irrep £1/2 nothing is so trivial since it is not 
possible to introduce a single group label for every com- 
ponent because of the bidimensionality of the irrep £1/2, 
however the general approach developed in the previous 
section will show its relevance in opening the way. 
One must take into account the fact that the V_ B (g) ma- 
trices are block-diagonal with dimensions 1,2,1. First 
from Eq. (|44l) we easily deduce directly that the first com- 
ponent of the two mutual partner eigenstates of the £1/2 

irrep (tp^ 1 / 2 (r), /i = 1,2), are mutual partner functions 
of the single group 2D irrep £ ! The same holds also for 
the last components ip 4 , \i — 1,2, but with a reversed 
association of the partners with the eigenstates. These re- 
sults can be relatively simply understood by considering 
the representation product '£3/2 © £1/2 = £, for i = 1,2. 
Much more intricate is the behaviour of the two cen- 
tral components of the two partner functions, i.e. 

VV, 2 /2 ( r )> M = 1,2 and (r), M = 1,2. Indeed, there 

one must resort to the full power of Eq. (|4"T)l to find 
the set of symmetrized envelope functions that form a 
basis for the 4D reducible representation £1/2 <8> £1/2- 
Using the products £i/ 2 ® £1/2 = A\ © A 2 © £ and 
E © £ = A\ © A 2 © £ it is possible to disentangle the 
central components into corresponding single group ir- 
reps. The full result reads as : 



±f 1/2 (r) = 



V2 



-0f(r) \ 
^i(r) + $f ( r )) 



V 



1 

1 

V2 



0F(r) 



W) 



frAi ( r ) 



*£(r)) 

iW) 



-2 



(53) 



/ 



where again a minus sign was added on <p A2 (r) . There 
are three different £ representations which appear in the 
expressions above. Let us identify clearly the different 
sets of B -functions with the help of the general theory : 



[Sl/2>-E 



J (r) 



(54) 



which helps to keep track of the different possible ori- 
gins of the £ representation. In total we see that with 
r = £1/2 we can associate the five independent UREF 
4> Al (r), 4> A2 (r)0f (r), ipf (r) and $f (r) (the partner func- 
tions (f>f (r) , ipf (r) , &f (r) are not independent, they are 
related to <f>f (r), (pf (r), $f (r) by a symmetry operation). 
Moreover multiplicities p» , r larger than 1 cannot be 



illustrated here because they would only occur in the cu- 
bic and icosahedral point groups which are the only point 
groups which are not simple-reducible [2l|. 

Despite of the fact that five independent UREF func- 
tions are now required for T — £1/2 (instead of the four 
independent envelope functions related to the original 
spinor ^ 1/2 (r)) we shall see in Section tVTl tlrat. thanks to 
the enabled SDR technique for spinors, the use of five in- 
dependent UREF's indeed leads to a maximum reduction 
of computing time and computing memory requirements. 



E. The symmetry of the Luttinger matrix elements 
in the OBB 



In this last subsection we shall shortly comment on 
the fact that simpler symmetry properties of the enve- 
lope functions goes hand-in-hand with simpler and more 
elegant expressions for the matrix elements of the k.p Ha- 
miltonian. It is possible to obtain scalar transformation 
laws for every matrix elements by considering again the 
expression of invariance of the ^-transformed Luttinger 
matrix under symmetry operations, and taking advan- 
tage of the explicit separation of symmetry operations 
into spatial and spinorial parts (c.f. Eq. (|2"5])). 
It can be shown that, in the OBB basis, the Luttin- 
ger Hamiltonian takes a simpler form, and also that 
every q, r, s operator can carry a single group irrep la- 
bel : p and q are scalar operators transforming under 
symmetry operations like A±, whilst the two operators 
(r, s) = (Of ,62) form a set of irreducible tensorial ope- 
rator (ITO) transforming like the partner functions of 
the irrep £, i.e. 



$(m ft E l9 (3D) 



M g 



(55) 



We refer the reader to l28| for more details. 



To shortly summarize the whole section [TVi let us em- 
phasize that the introduction of the concept Optimally 
symmetrized Bloch function Basis (OBB) corresponding 
to the basis reducing the representation of symmetry ope- 
rations in spinorial space allows to systematically decom- 
pose envelope functions into a minimal set of minimally 
coupled UREFs (Ultimately Reduced Envelope Func- 
tions) carrying only single group irrep labels which spe- 
cify their individual transformation properties. It should 
be pointed out that the coupling in question is through 
symmetry operations, the Hamiltonian will still couple 
these functions together, following the single group irrep 
multiplication table. 

From now on we shall leave out the underscore nota- 
tion on matrices and spinors, assuming that the reader 
is familiar with the implicit nature of these basic objects 
of the theory. 
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V. SPATIAL DOMAIN REDUCTION 
TECHNIQUE FOR SCALAR FUNCTIONS 

In this section we develop a technique which we call 
spatial domain reduction (SDR), and which allows one 
to maximally reduce the geometrical domain of solution 
of eigenproblems linked with partial differential opera- 
tors for scalar functions. It is the second essential ingre- 
dient of the MSRF formalism. The technique is closely 
related to the general theory developed in ref. [42j], but 
identifies from the outset the essential elements of the 
symmetric domain and the separate role of each special 
boundary and singular point. Other elements of the un- 
derlying theory, based on symmetry adapted basis func- 
tions, can be found in (43|. Although this SDR method is 
applicable to any point-group symmetry, we shall again 
confine ourselves here on the Czv case, keeping a rather 
intuitive and pedestrian approach, sufficient for grasping 
easily the method, and using it. A full fleshed mathemati- 
cal treatment encompassing these developments has also 
been set up, but will be given elsewhere. In the following 
section (Sec. lVI|) we shall also show how to go beyond the 
scalar case and reduce the spinorial case, but this requires 
a simultaneous use of OBB, consistent sets of UREF's, 
and the SDR. The formulation will apply equally well to 
a spatial 2D problem or 3D problem, i.e. a quantum wire 
or a quantum dot (c.f. Eq. ([To]) ). 

An important point to stress is that the SDR technique 
is general enough to simplify the treatment in the same 
way for all real space methods provided that the discreti- 
zation scheme respects symmetry. We have actually used 
a FE method, but the procedure would be left unaltered 
for FD methods : one only needs to identify to which 
domain a given node belongs to separate the blocks in 
Eqs. ([551) and ([751 to follow. A completely similar ap- 
proach would even hold in the case of tight-binding (TB) 
methods provided one is able to identify the spatial lo- 
calization related to each element of the basis. Further 
along such a line of thought other versions of the SDR, 
OBB, and MSRF techniques could be developed for Fou- 
rier space formulations, using analogous separations of 
the Fourier part and spinorial part using a symmetrized 
basis in both space (an example of symmetrized plane 
wave superpositions covering all sectors of the Fourier 
space is given in 

Let us now consider any arbitrary scalar Hamiltonian 
H displaying the C^ v group symmetry and related to a 
2D structure like presented in Fig. [7J We shall detail in 
the next three subsections the three different steps of 
the method which are : 1) decomposition of the spatial 
domain into minimal disjoint sub-domains, 2) identifica- 
tion of the minimal set of independent sub-domains for 
every irrep (the reduced domains), 3) computation of the 
correspondingly reduced Hamiltonian matrices on the re- 
duced domains. The procedure is systematic, and allows 
to treat non-trivial cases like the minimal domain for the 
2D irrep E in G$ v . 



A. Decomposition in disjoint sub-domains, and 
representation of symmetry transformations by 
permutation matrices 

We consider a computational domain T> with C^ v sym- 
metry (Fig. [5]) , and assume that eigenstates ip are repre- 
sented by a vector of dimension N, where each vector 
component ipi represent the value of the wavefunction 
ip{vi) on every mesh node i at position in a FE or 
FD approach. In a TB approach ip^ would represent the 
weight of an orbital at each site. 

Since the domain T> and the system Hamiltonian are 
C^-symmetric, we can set without restriction a symme- 
tric spatial basis with respect to every operation g of 
Cs v . From a real space numerical point of view, it im- 
plies the use of a symmetric FE or FD mesh. Let us 
then decompose the full domain T> in thirteen disjoint 
parts (sub-domains) by separating the six interior parts 
£j, i = 1, . . . , 6 in-between the three symmetry axes (c.f. 
Fig- HI j the six borders Bi,i = 1, . . . , 6 between the inter- 
ior parts, and the remaining central point C (this domain 
decomposition must be performed in such a way that the 
domains are always the image of one another under every 
symmetry operations). 

Let us now formulate the discretized form of Eq. (fT3| un- 
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Figure 9: Decomposition of the spatial domain in disjoint 
parts 

der every change in coordinates corresponding to a sym- 
metry operation g £ Cs v : 

r- = O^r* = r i=7rs(i) (56) 

This condition obviously defines a set of permutations 
■Kg of the basis nodes which do form a representation of 
the symmetry group C^ v . Let us now assume that within 
the thirteen sub-domains identified previously the node 
numbering is always a perfect image of the other similar 
subdomains under any symmetry operation (this does 
not restrict generality : a further inner-domain permuta- 
tion could be added to eventually describe arbitrary non- 
symmetric node- numbering between the sub-domains). 
Under this simplifying assumption all the permutation 
matrices become a block matrix permuting only the sub- 
domains, and leaving unchanged the internal subdomain 
node numbering. We are thus led to consider 13 x 13 "do- 
main permutation matrices" V g , with the identity within 
each block, and defined as Vij(g) = 8^ To give a 
concrete example let us represent the matrix V Gvl corres- 
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ponding to the vertical mirror operation a v \ (see Fig- [7j) , 
as well as its transposed action on the set of subdomains : 
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where the big matrix is All the other operations 

and associated matrices can be straightforwardly written 
down in an analogous way. Considered as simple 13 x 13 
"domain permutation matrices" , the set of V g permuta- 
tion matrices obviously forms a reducible representation 
of the symmetry group C^v ■ 

The next step is to decompose an arbitrary wavefunc- 
tion ip( r i) hrto blocks corresponding to the domain de- 
composition : 

/ V>e \ 
/ Vbj \ 

4>s 1 
V>s 2 
Vs 2 

V>B 3 

tp = Vs 3 (58) 
Vs 4 

V>B 5 
V>S 5 
I V>B 6 I 

\i>s e I 



where ipc is the central value, ipBi the subvector of 
values along the border B\ (Fig. [9]), Vsi the subvector 
representing the values of the wavefunction inside the 
sub-domain Si, etc... In the numerics it is sufficient to 
group the nodes by domains. Note that in this approach 
we do not need to refer to the dimensionality of the 
structure, i.e. the technique will work equally well for a 
2D spatial problem (quantum wire), or a 3D problem (a 
quantum dot). 



It is also important to note that, according to Eq. (fT5[) . 
a passive change in coordinates described by g € C% v is 
now equivalent to perform a V g permutation of the 13 
subparts of the wavefunction : 



(59) 



The decomposition (|58|) of the subparts of the wave- 
function ip, and the identification of its transformation 
properties (|59p . could seem simplistic, but they are in 
fact an essential step for the formulation of the theory, 
allowing in the foregoing section the identification of the 
independent geometrical parts of the wavefunctions, as 
well as, much later, a compact solution of the problem 
in terms of reduced geometrical variables, i.e. on a mini- 
mal set of independent sub-domains for each irrep of C^ v . 



B. Minimal set of independent sub-domains 



Let us construct the projection operator on the 
partner function /i of a given irrep T of C% v , operating 
on wavefunctions. Group theory provides a closed form 
expression for where we can use the set of V g permu- 



(57) tation matrices just defined 



^frE(^)); 

11 geg 



(60) 



Therefore any arbitrary function on T>, represented by a 
vector ip, can be written uniquely as a sum of symmetri- 
zed components ip^ : 



(61) 



In particular, a function ipj^ transforming like the partner 
function ji of the irrep T, is necessarily an eigenstate of 



the projector with unit eigenvalue : 



(62) 



This condition allows us to identify for every irrep the ne- 
cessary relations between the values of ip^ on the different 
subdomains of T>. For example, for r = A\ or T = A2, 
the condition (|6"2")l leads to 
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(63) 



where the superscript T = A\ (respectively A2) is impli- 
cit on every components tpc > 4>Bi > " • of ip Al (respectively 
Tp A2 ). The above results (|63[) for A\ or Ai are very inter- 
esting and deserve two comments : first they show that 
each wavefunction is characterized by its values on a mi- 
nimal subset of subdomains (which is not the same for 
all irreps), second for A\ and A2 it should be possible 
to characterize the wavefunction essentially on one sixth 
of T>. Note also that the borders are not necessary for 
A2 since Eq. (1551 shows that ip A2 must necessarily be 
zero on the symmetry axes between the six main subdo- 
mains ! For numerical efficiency one is therefore lead to 
seek a solution for each irrep separately, on the minimal 
set of subdomains only, and this will be considered in 
detail in the next subsection, after the treatment of the 
less trivial case of the 2D irrep E. 

A novel aspect of the 2D irrep E is that there are a 
few equivalent matrix representation of Cs v , and a spe- 
cific choice must be made beforehand. Since with our 
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axis choice the 3D rotation matrix 5ft do already form 
a reduced 3D representation of the C$ v group, it is al- 
ready in the desired block- form (A\ © E). For simplicity 
we shall take this specific matrix representation for E 
(explicitly presented in the appendix [Al . One advantage 
is that this representation is real, another one is that 
the matrix representing a v \ is diagonal, hence the first 
and second partner functions will be respectively even 
and odd with respect to <j v \. The last advantage is that 
conditions of parity with respect to u v \ correspond to 
local Dirichlet or Neumann boundary conditions (if ins- 
tead one would have chosen to diagonalize the rotations 
C 3 one would have got complex wavefunctions with non- 
local 27r/3-periodic boundary conditions as a function of 
the azimuthal angle). 

For each partner function one can now build a different 
projection operator, called and P® according to 
Eq. (|60[) . We can however work with only one partner 
function, because the independent parameters are rela- 
ted to the irrep and allow to re-built the second partner 
function from the first. For each partner function one 
finds : 



This leads to the following full expressions for the first 
partner function as a function of the values of the two 
partner functions on the first sixth of the domain : 



o 

01, Bi 

-2^1, 8 2 
-^l.Si-^^a.Si 



o 
o 

-\/3 01,B 2 
-|V-2, Sl + ^^l, Sl 

o 



v"3 



■ -01,81 



. s/3 



■01,51 



V 



-V3 0>1,B 2 



J 

(66) 

For the 2D irrep E, rjjf and ipf are obviously less in- 
tuitive than ip Al and %p A2 . However the important point 
is that we have found a systematic procedure to derive 
them which even allows to find straightforwardly novel 
analytical results. For instance we can see in Eqs. 
and dHJ) that 



V 





^l.Si 
01, Si 
■01, B 2 
01, S 2 
-|-01,Bi 
V'l.Si — V , 1,S 2 

-2^1, B 2 
■01.Si-01.S2 
-^l.Bi 

0i, s 2 

■01, B 2 

0i, Si 



and 



n 


i>2,S 1 
-02, B 2 
V>2,S 2 
02, B 3 
-4'2,S 1 +02. S 2 


02, Si -02, S 2 
- -02,63 
—02, S 2 
— 02, B 2 
— 02,5i 



/ 

(64) 

We see that the corresponding reduced domain is now 
larger than a sixth, and its surface/volume exactly cor- 
responds to a third of the full domain. It is seemingly 
not possible to rebuilt the full function from a smaller 
domain as for the non-degenerate Ai irreps. It is however 
possible to bypass this limitation, and to still reduce to 
one sixth of the full domain, by noting that, since the 
first and second partner functions are not independent, 
one could retain simultaneously the two of them, but only 
on the first sixth S\. This amounts to replace the "01,52 
variable of the first partner function by the -02, Si variable 
of the second partner function, related to the first inter- 
nal domain ("first sixth"). It is possible to relate these 
variables using the a v s symmetry operation. Using the 
transformation rule (1301) we find : 



01,S 2 -> [D E (<T v3 )] n 0i,5i + [ DE {< j v3)] 21 02,5i 

= -^l,Si + ^02,5i 

02, B 2 -> V30i,e 2 
V3 . 

1p2,B 3 -> — 01,Bi 

1 / , V 5 , 

02,S 2 "J" 2^ 2 ^i + —VlSi 



(65) 



Iplfi = 02,C = 



(67) 



This is non-trivial and had not be noticed before, even 
with correct numerical results at hand. 



Reduced domain and non-trivial boundary 
conditions 



We are now going to use the identification of a mi- 
nimal number of independent part in the wavefunction 
to reduce the size of the problem at hand. In addition 
we show how non-trivial boundary conditions naturally 
occur at the border of the reduced domain. 

For every irreps we may collect the minimal set of in- 
dependent parameters and define a reduced vector ip^Red 
on the corresponding reduced domain. For the irreps Ai 
and Ai they are given by : 




and for the irrep E by 



^Red 



and 



VRed 



Vte, (68) 




(69) 



A compact way to summarize Eqs. (f6"8"|) - ([6"9"]) is to in- 
troduce for every irrep F a set of rectangular reduction 
matrices , fi — 1, . . . dr such that the reduced vectors 



can be written 



Red 



(70) 
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The reduction matrices S£ will be used in Section fVE II 
to obtain the reduced Hamiltonian corresponding to the 
independent variables. Note that the Sj^ matrices, obtai- 
ned here in a rather pedestrian way, are not unitary (for 
another approach see e.g. ref. [45| and references therein). 

Let us now shortly look at the natural boundary condi- 
tions which occur at the edges B\ and Bi of the reduced 
domain on one sixth S\. For the ID irrep A2 the boun- 
dary conditions are obvious from Eq. (155)) and are of the 
Dirichlet type, since the wavefunction naturally vanishes 
on the border : 

< 2 - < 2 = ° (71) 

For the ID irrep A\ it is not very difficult to show that 
the corresponding boundary condition is of the Neumann 
type, and that the normal derivative vanish on the bor- 
der. To show this we note that the behaviour of the 
normal derivative on an edge Bj can be connected to 
the difference between the subdomains on each side : 
~ (ipB^-i ~ V'g 1 ) where we assume the definition 
V'jBo = ^bI- Therefore the normal derivatives ip'^ 1 and 
^'b2 van ish, i- e - 

^=^=0 (72) 

These intuitive results for the ID irreps A\ and A2 were 
clearly apparent in the numerical behaviour, have been 
easy to enforce numerically when computing on the redu- 
ced domain. Furthermore they correspond to local boun- 
dary conditions in a real space formulation (FE or FD 
methods), hence they would not destroy the sparsity of 
the Hamiltonian and mass matrices. Finally to alleviate 
the need for node re-numbering the so-called "penalty 
method" can be used to enforce the Dirichlet boundary 
condition 46]. 

Let us now look at the natural boundary conditions 
for the two partner functions linked with the 2D irrep 
E. We shall see that these natural boundary conditions 
become not so trivial, and could hardly be guessed by 
just looking at the numerical result. It may seem easy to 
directly read valid boundary conditions from Eqs. (|64|) 
or dSBJ), i.e. V'i £ Bl = i>' x % = and = ^e 4 = 0, ho- 

wever these boundary conditions are not on the reduced 
minimal domain, and using them we would only take ad- 
vantage of C s symmetry only, not the full Cs v symmetry 
of the problem. 

Therefore let us look more closely at what happens 
at the border of the reduced domain in the first partner 
function in Eqs. or (|B"6"j) . Clearly on the border B3 
we have 

^B, = -^iV (73) 

This beautiful analytical result is not very convenient to 
use : from a numerical point of view it does correspond 
to a non-local boundary condition (coupling far apart 
points), and a pedestrian implementation would destroy 



the matrix sparsity pattern linked with nearest neighbor 
mesh nodes. 

For the second partner function a similar path allows to 
obtain the corresponding non-local boundary condition : 
this time the normal derivative of the function with res- 
pect to the boundary will be involved. Since along edge 
Bj the normal derivative of the second partner function 
i/> 2 E b i s related to the difference of the function across the 
boundary, V'^B- ~ W^b-_i ~^2Bi)' ^ is easy to see from 
the second equation in Eqs. (IMl) or (|BB|) that the function 
if>2 necessarily obeys the following boundary condition on 
the border B3 : 

The non-trivial character of this boundary condition is 
further evidenced by the fact that on the border B3 the 
value of the function ?/>f' is not vanishing, by contrast 
with the value on the border B\. 

The problems linked with boundaries have been stu- 
died in detail by ref. [42[ with the help of a general ap- 
proach, valid for any symmetry group, but also applied 
specifically to C-$ v as an example. Eqs. (|72p and (J7TJ) were 
straightforwardly found, but not directly Eqs. (l67j) . (f73|) - 
(|74[). because their approach considers only coupling the 
two partner functions on the boundaries of the so-called 
"fundamental symmetry cell" <Si, as could be read di- 
rectly from Eq. ([66]) . 

To summarize, in the present section we have found 
the boundary conditions (|71ti74l) which do allow to solve 
numerically in an unambiguous way the eigenproblems 
related to each irrep on a minimal domain whose size vary 
from one sixth to one third. However for the 2D E repre- 
sentation they are both not very trivial and non-local in 
space. Later, in subsection IV El we shall see that the use 
of SDR-reduced Hamiltonians will avoid the explicit use 
of boundary conditions whilst keeping the natural spar- 
sity patterns, i.e. it will be possible to solve all issues at 
once. 

D. The structure of the full Hamiltonian 

Having developed the SDR technique for every type of 
wavefunction, we shall now consider the form of the full 
Hamiltonian (i.e the Hamiltonian on the full domain). 
This will enable us to reduce the full eigenvalue problem 
to distinct reduced eigenvalue problems for each type of 
irreducible representation. We recall that we deal for the 
moment with a scalar Hamiltonian. In the case of C% v 
such a Hamiltonian can be seen as an invariant 13 x 13 
block operator, operating on the decomposition (|58p of 
the subparts of the wavefunction ip. Then using group 
theory it is easy to construct the most general 13 x 13 
block-form matrix which would satisfy the following two 
conditions : 1) respect of the sub-domain connectivity 
induced by differential operators, and 2) invariance with 
respect to the symmetry operations of the group (i.e. the 
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full Hamiltonian operator belong to the irrep A\). The 
discrctized form of the Hamiltonian will be represented 
by a similar block-matrix. 

The first condition, concerning the most general Ha- 
miltonian on the full domain respecting the connectivity 
between sub-domains, can be easily enforced : the Hamil- 
tonian must have a block-form corresponding to the the 
most general 13 x 13 matrix H c respecting the connecti- 
vity, which is written down by crossing out off- diagonal 
elements that would couple non- contiguous domains in 
real space. Here, since we use only first order FE elements, 
we can also take into account the further simplification 



that the center C couples with every edge Bi but not with 
the internal domains <S, (but this restriction can be easily 
lifted). 

The second condition, concerning the invariance, can 
be satisfied by enforcing the fact that H must be inva- 
riant under projection on the irrep A\ (the projector on 
A\ for operators^ a "superprojector" , is built analogously 
to (|60p , where the linear operator V g is replaced by a li- 
near "superoperator" , whose action on an operator O is 

One can then find that the corresponding scalar Ha- 
miltonian matrix H is necessarily of the form 



H = H? 1 = 



He H c .b 1 -He, Si Hc,B 2 



Hc,s% Hc : B t H c ,sx Hc,B 2 Hc.Sx \ 
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H Sl 
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B\ ,B 2 


S lt B 2 
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V Ht, Sl 


Sl.Si 
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H S t ,B 2 


Hs 



In this equation the square diagonal blocks Ha are simply 
noted Hi and represent the restriction of the Hamilto- 
nian operator within every subdomain. The off-diagonal 
rectangular blocks Hij = fft couple the i and j sub- 
domains. One can see in Eq. (|75[) that there are many 
redundant blocks : indeed this is a consequence of symme- 
try and it was one of our primary goals to identify them. 
FD or FE schemes will automatically generate such a 
block structure for the Hamiltonian, provided the mesh 
does respect C^ v symmetry. Each block can be either 
computed numerically only once (to save computer time, 
but it requires some care), or can be obtained by sam- 
pling a posteriori the blocks in the full discretized Hamil- 
tonian. To this end it is sufficient to tag every node with 
its domain. 

It should be pointed out that the reduction technique 
presented in the next section allows to go even further in 
the reduction by building the smallest possible reduced 
Hamiltonians, where every Hamiltonian block appears a 
minimum number of times. The proposed procedure will 
amount to analytically pre-diagonalize by block our full 
Hamiltonian using group theory. Then the final numerical 
diagonalization of each block on the diagonal will directly 
produce the independent relevant parts of each type of 
eigenfunction without redundancy. 



E. The reduced Hamiltonians and explicit results 
for the C'sv group 

For every irrep, it is useful to obtain a corresponding 
reduced Hamiltonians related to the minimum number 
of independent part on the sub-domains. To this end the 
rectangular reduction matrices <S^, fi — 1, . . . dp defined 
earlier in Eq. ([70)) allow one to write : 

H r Red = ^~ 0s£) f H Sl (76) 

where du e d represent a "normalization" factor, chosen 
conveniently so that the diagonal blocks Hg 1 of the Ha- 
miltonian, corresponding to the interior parts Si — S§, 
would be unaffected by reduction. By definition the num- 
ber of blocks in such a reduced Hamiltonian is minimal 
and its solutions correspond to a given irrep T of the 
symmetry group of the full Hamiltonian. It should be 
pointed out that it is enough to solve only once within 
each irrep T to obtain all eigenvalues and reconstruct all 
eigenspaces by using symmetry transformations. This is 
also reflected in the fact that -ffjj ed is independent on fi. 

In a FE formulation of the eigenvalue problem one 
should also take into account that the original differen- 
tial equation in real space is mapped to a generalized 
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eigenvalue problem 



(77) 



where H and M are respectively the so-called "stiffness" 
and "mass matrices", and E^,n = 1, ... are different ei- 
genvalues all labelled by the T irrep. The stiffness ma- 
trix H correspond to the Hamiltonian expressed in the 
FE basis and has also the block- form presented in (|75|). 
Similarly the mass matrix M which represent the non- 
orthogonality of FE basis is also invariant (i.e. is associa- 
ted with the A\ irrep) and has the same block-structure 
as H . The scalar product in the FE approach simply re- 
duces to {ip\(j>) = i\) + M (j), where -0 and <f> are the vectors 
corresponding to the coefficients of the decomposition on 
the FE basis. Hence in the FE approach one also needs 
to define a reduced mass matrix similar to (1761) : 



M Red 



1 



dRe 



so that the set of reduced problems then reads 

HRed ^Red 



E T M T Red ip T Red 



(78) 



(79) 



where the factor dn e d cancels. It may however be inter- 
preted : since both matrices M and M Red can be conside- 
red as measures involved in the definition of scalar pro- 
ducts of vectors of type ip and ipR e d appearing in (|77|) 
and (|79j) respectively, the value of dn e d corresponds to 
the ratio of the measures of the full and the reduced do- 
main (but would be unity for unitary Sj^ matrices, see e.g. 
ref. [4a | and references therein). In our C% v case our pro- 
cedure led naturally to the value du e d = 6, which indeed 
is linked with a domain reduction by a factor 6. Eq. (|76[) 
shows that there are as many separate eigenproblems as 
the number of irreps in the group Q . In a FD scheme the 
mass matrix is reduces simply to the identity. 

We are now going to review the specific forms of the re- 
duced Hamiltonians for our C$ v specific example. Similar 
forms will also hold for mass matrices. 



Reduced Hamiltonian for A\ 



For A\ the reduced Hamiltonian reads 



H 



Red 




(80) 



The factor 1/2 appearing in front of the Hb 1 and H& 2 in 
the diagonal can be interpreted simply in a FE scheme : it 
is equivalent to fill the matrix with only the contribution 
of the FEs on the interior side between the symmetry 
planes a v \ and a V 2- It can be proved that such a pro- 
cedure implicitely enforces a Neumann condition (with 
a zero normal derivative on the interface), evidencing 
clearly that indeed the boundary conditions correspon- 
ding to A\ symmetry are implicitely incorporated ! The 



same intuitive argument can explain the factor 1 /6 for the 
central block He - In Fig. [TO] we show the ground electro- 
nic eigenstate with A\ symmetry in C-^v quantum wire. 
The highlighted sixth correspond to the reduced spatial 
domain used in the numerical solution. 




Figure 10: Eigenstate with A\ symmetry 



2. Reduced Hamiltonian for A2 

For the A2 irrep, the reduced Hamiltonian will ob- 
viously have only one block (see Eq. ([581 ). 



H Red - H Sl 



(81) 



The eigenstates then naturally vanish on every border 
(Dirichlet condition) and again one only needs to solve 
on the first internal sub-domain. Fig. [TT] displays such 
an eigenfunction, found as a highly excited state in C% v 
quantum wire. For convenience in our numerics we still 




Figure 11: Eigenstate with A2 symmetry 

use the same grid as for A\ and treat the border nodes, 
but "cross-out" the border points and enforce the Di- 
richlet boundary condition with the "penalty method" . 
This avoids a cumbersome node-renumbering task bet- 
ween the two irreps A\ and A2 , at a negligible numerical 
cost. 



3. Reduced Hamiltonian for E 

Finally, let us consider the more complicated degene- 
rate irrep E and apply the same procedure. We obtain 
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the following reduced Hamiltonian 



VI. THE CASE OF SPINORIAL SETS OF 
FUNCTIONS 



n Red 



H Sl H Sl .e 2 
C.C. 2ff B . 



(82) 



corresponding to the "mixed" reduced partner wavefunc- 
tion (l69l) . The s/3Hs lt B 2 block represent the non-trivial 
coupling between the two partner functions (coupling of 
"01. Bi with -02, Si)- As explained before, in the same way 
we could have started with the second partner function 
and obtain a similar reduced Hamiltonian with respect 
to a different but similar set of reduced variables (see the 
dependence of on fi in Eq. ([701) 1. but we do not need to 
do so because the new Hamiltonian would carry exactly 
the same information, and starting from the reduced va- 
riables in i/jft ed one is already able to reconstruct the full 
domain vectors V'f and ipf given by Eqs. (|66| or Eq. (|64|) 
and its corresponding form for ip^ ■ As an illustration we 





(b) 



Figure 12: Eigenstate with E symmetry ((a) partner func- 
tion 1 (even), (b) partner function 2 (odd)) 



show in Fig. [T^] the corresponding degenerate even/odd 
partner eigenfunctions linked with the first excited elec- 
tronic state in a C 3v quantum wire. 

From the numerical point of view the advantage of 
the "mixed" reduced partner wavefunction (|69l) on one 
sixth of the domain is that one can use as elementary 
variables at every node the two values corresponding the 
first and second partner functions. This allows to keep 
the same sparsity pattern linked with the connectivity 
on one sixth, but with blocks of two variables. In this 
way node-renumbering can be avoided also for the E re- 
presentation. 

Among the important points of this section, let us men- 
tion that now there is no need to incorporate explicitely 
any of the non-trivial boundary conditions in the set of 
reduced problems (|79)l with H^ ed (and eventually Mjj ed ) 



given by Eqs. (JHOJ) or (JHU) or flS2J|, they will be auto- 
matically satisfied, whilst the hermiticity of the matrix 
problem will be preserved. 



In the last sections, we presented two essential me- 
thods allowing to choose an optimal bloch function basis 
for a spin dependent problem (Sec. IIV|) and to perform a 
spatial domain reduction for a spinless problem (Sec. |Vj 
. A combination of these two approaches allows to ea- 
sily reduce spin-dependent problems, as we show in the 
present section. For convenience we shall restrict again 
our discussion to the C 3v case, although the procedure is 
general and quite simple : one starts with the spinorial 
eigenstates (|5"Tj) - (f5"3")) . labelled by the double group irreps 
r, and, using the transformation defined by Eq. (|47p . 
one can construct the corresponding vector of UREF's 
for each irrep T = 2 E 3 / 2 , 1 E 3 /2,E 1 / 2 - Since every UREF 
can be considered as a scalar function labelled by a single 
group composite label (irrep + partner labels) , one is able 
to treat them with the SDR technique by dividing each 
one into 13 pieces corresponding to the 13 subdomains. 
The reduction must then be carried out simultaneously in 
the sets of UREF's linked with every irrep T, by selecting 
the minimum number of independent parts. 

This composite procedure leads to the following three 
reduced sets of UREF's : 



^3/2 



Red 



'Red 
''Red 
'Red 



( 



^3/2 

Red 



'Red 
''Red 
'Red 



(83) 



VRed 



(84) 



<t>led 

mM 

VRed 
^Red 
^Red 
V Vied J 



To simplify the notations we have omitted in Eq. (|83[) 
the implicit double group label 2 E 3 / 2 which should also be 
beared by every component function 



'Red' 



Red and 

E, 



3/2 

Red 



'Red 

and 



of the reduced spinor ^Red, 2 ( an d similarly for %jj 

^Red)- Despite the fact that in Eq. (|84| the degenerate 
irrep E appears three times in the UREF decomposition 
(c.f. Eq. (f53"|)). we outline that the corresponding UREFs 
't'Red' ®r edi and Plied &Te distinct functions and corres- 
pond to independent variables. Using the reduced spi- 
nors (|53")1 and ((8~4")) one can construct in each case the cor- 
responding Sf matrices (analogously to Eq. (1701 ) which 
allow to reduce as in ((75)) the full Luttinger Hamiltonian 
on a sixth of the structure for each double group irrep T. 

The last step to achieve the full reduction of the Ha- 
miltonian on one sixth of the domain is to construct the 
most general 52 x 52 block-form valence-band Luttinger 
Hamiltonian operating on every subdomain separately. 
We shall take advantage of the fact that, as discussed in 
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Sec. IIV1 the p, q, r, s operators can be decomposed with 
respect to the single group irreps in a similar way as for 
the spinorial eigenstates. In particular, for C$ v , p,q are 
Ai-type operators whilst the couple (r, s) forms an ITO 
transforming like the E irrep. 

Hence since p and q are Ai, their most general form 
is identical to the form of the scalar Hamiltonian ap- 
pearing in Eq. (|75p . whilst the block form of r and 
s must be specifically computed using the restriction 
that they are partner operators in E. To this end one 
must simply restart from the 13 x 13 connectivity block 
matrix H c and use the projectors on E : H^ ,z — 
T§T £ 9 es Af*G?) 'PgHcVg 1 , with i = 1,2 for r, a respec- 
tively. In this way the full 52 x 52 block Hamiltonian can 
be constructed. 

The reduced Hamiltonian and mass matrices H^ ed and 

Mft ed for each irrep T = 2 E 3 / 2 , l E 3 / 2 , E 1 / 2 can be obtai- 
ned using the new reduction matrices Sf resulting from 
Eqs. (|8"3"f and The form of these reduced Hamilto- 
nian is discussed in the foregoing subsections. 



C 



-s Sl 
rs, 



.V3r 



and depend only on the r and s operators. These terms 
are a consequence of band coupling between different en- 
velope functions which is a feature of the Luttinger Ha- 
miltonian (j4]) . This band coupling is fully compatible with 
the wavefunction decoupling during symmetry operations 
achieved by the OBB (c.f. Eq. If one would neglect 

band coupling, i.e. band-mixing, the reduced Hamilto- 
nian would be a block-diagonal set of decoupled scalar 
Hamiltonians. 



Since r and s are not self-adjoint operators, we have 
Cai-e 7^ Ce-a > bvrt the natural correspondence is 



Ca z -i 



where tv, and Si 



are the blocs appearing in the r, s operators (c.f. Eqs. (JS7J- 



A. Reduced Hamiltonian for the non-degenerate 
irreps l E 3/2 

For the non-degenerate irreps l E^/ 2 i i = 1,2, the redu- 
ced Hamiltonians read 



H E ' 3/2 

n Red 



TT E 3/2 

n Red 



HrUp + v) C A ^ E (r lS ) 







2Hi ed {p- 



h.c. 



q) CE-Ai(r,s) 



(85) 



H Red(P + <l) C Al - E {r,8) 







h.C. H^ d {p + q) 



(86) 

On the diagonal we find that H^ d , H^ d , H^ ed have a 
form exactly similar to the reduced Hamiltonian for the 
spinless conduction band problem (|80I81|) and (|82|) . and 
depend only on the p + q and p — q differential operators. 

The non-diagonal coupling terms Cr 1 -r 2 , — ^2 
or E have the following forms 

C E -Ai (r, s) = 

\ r B 1 ,C |r Bl r Sl,Si r B l ,B 2 

r Si,Bi r Sl r S 1 .B 2 

s Si,i3i s Si s Si,B 2 

2l'B 2 ,C r B 2 .B 1 +V3 SB 2 ,B 1 '"B 2 ,S 1 +V / 3 SB 2 ,Si 2 r B 2 

(87) 



It is apparent in the structure of the reduced Hamil- 
tonians (1551 and ([5d| that they are linked to each other 
by time-reversal invariance. Nevertheless in the case of 
quantum wires one should pay attention to the fact that 
the parallel wavevector k is reversed by time reversal, so 
that for the same k- value the Hamiltonians (|85[1 and (|55|) 
give rise to different eigenstates. In quantum dots the ei- 
genstates of ([55)1 and ([55)1 form Kramers degenerate pairs 
and it is enough to solve for one of them only. 



B. Reduced Hamiltonian for the degenerate irrep 

El/2 



Let us now discuss the most complex case, the 2D E\/ 2 
irrep, with the same reduction technique. We obtain the 
reduced Hamiltonian : 
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n Red — 2 



2Hf, 



Ap-i) 



, \ 

-^C E -E 



V 



775 ^ 



(89) 



It is interesting to note that the coupling terms Ce-Ai and Cj^—Eii = 1,2, are the same operators which appeared 
in the previously reduced Hamiltonians (pj|) and (|86p . However there is only one additional block, Ce-e, given by 



C E -e{t, s) 



2 r Bl 
r Sl,B 1 



''Si 
-SSi 



^ '•B2,B 1 +V / 3sb 2 .Bi r B 2 ,Si-V / 3sB 2iSl 
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C. MSRF post-symmetrization technique 

In the previous subsections we have simultaneously 
applied the OBB and the SDR techniques, what we 
call the MSRF technique, to decompose totally the 
envelope functions into "Ultimately Reduced Envelope 
Functions" (UREF) on a minimal domain, for a problem 
involving spinors and subject to double group symmetry. 
As result of the explicit separation of the spinorial 
and spatial part, in the double group Hamiltonians, 
all UREF's were shown to depend only on the single 
group irreps, for C^ v the A\, A2 and E irreps. We have 
also obtained reduced Hamiltonians and showed that 
they have in this formalism fully symmetrized matrix 
elements. The reduced Hamiltonians are block-diagonal 
and their structure is as follows : diagonal blocks 
correspond to scalar Hamiltonians H T (p ± q) and off- 
diagonal coupling blocks Cr 1 ,r 2 (r, s) which only depend 
on the operators (r, s) and where the irreps Ti and 
T2 are single group representations which correspond 
to a coupling between UREF's with symmetry Ti and T2- 

One may however wonder about the complexity of the 
MSRF technique, which could be cumbersome numeri- 
cally because of the need to recode existing k ■ p codes. 
Actually it is not absolutely necessary to modify an exis- 
ting k ■ p code to benefit from most of the advantages of 
MSRF. We shall therefore conclude this section by out- 
lining how the MSRF symmetrization technique can be 
used as a "data post-processing technique" . 

The steps of MSRF post-symmetrization can be des- 
cribed are as follows : 

1. Identification of the main symmetry elements : a) 
the symmetry group by taking the common sym- 
metry elements between the k ■ p Hamiltonian and 
the heterostructure, and b) the subdomains for the 
SDR decomposition and buildup V g matrices. 

2. Construction of the OBB, this requires the follo- 
wing steps : a) identification of the set of V_(g) ma- 
trices linked with the original k-p envelope function 
formulation and check of Eq. (|2"9"|) . b) analytical 



block-diagonalization of the representation V_(g), 
and c) definition of the OBB via Eqs. (j34l [37l I38 |) . 

3. Numerical computation of all eigenstates of inter- 
est with existing k ■ p code and interpolation of the 
envelope wavefunctions on a symmetrized grid cor- 
responding to the SDR subdomain decomposition 

4. Symmetry classification of numerical eigenstates. 
This requires the following steps : a) for all ei- 
genstates perform numerically the change in ba- 
sis towards the OBB, b) symmetry classification of 
the interpolated eigenstates by computing the ex- 
pectation value of the projectors on every irrep T 
(identify also accidental degeneracies), c) identify 
partner functions belonging to the same eigenspace 
with the help of symmetry operations (and ave- 
rage the corresponding eigenvalues if there would 
be a slight numerical lifting of the degeneracy), d) 
change in basis within every eigenspace to further 
fix the form of the representation matrices linked 
with partner functions, e) eventual check the trans- 
formation laws P0|) . 

5. Construction of the UREF's using Eq. (gTJ), and 
eventual check of UREF's transformation laws (|49l) . 

Note that the order of some of these steps can be inter- 
changed. All analytical steps can be easily computed and 
automated on a personal computer. In [45l 1471 ] we have 
put into practice the automatic recognition of mode sym- 
metries, as well as the techniques of postsymmetrization 
and postconstruction of UREFs (called URCFs in the 
case of [iBj). 

We see here that MRSF symmetrization can be viewed 
as a "data post-processing technique" , but which still 
does allow to benefit from all subsequent advantages. 

In the next section, we shall demonstrate in practice 
the potential of MSRF by computing matrix elements 
linked with optical transitions. We shall show that 
one can obtain in our C^-example novel analytical 
expressions for optical polarization anisotropy which 
were far from obvious beforehand, and which go beyond 
the standard use of selection rules for eigenstates. 
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Typically such specific advantage of MSRF is due to the 
formulation in terms of UREF's, and the possibility to 
take advantage of selection rules at the intermediate level 
of UREF's. Another big numerical advantage of MSRF 
symmetrization will be that in any type of subsequent 
matrix clement computations the reduced domain can 
also be used, speeding tremendously computation. This 
can be especially important in heavy cases, for example 
when computing multidimensional Coulomb integrals 
for exciton complexes in quantum dots, or any complex 
object like polaron, etc. . . Like any symmetrized theory, 
MSRF enables the use of minimal coupling scenarios on 
a wide scale, even for possible fine-structure splittings 
due to symmetry breaking effects of all kinds. 



VII. SELECTION RULES AND MATRIX 
ELEMENTS WITH THE MSRF FORMALISM 

In this section, we apply the MSRF formalism to 
computation of matrix elements of operators. Wigner- 
Eckart Theorem (WET) allows to obtain selection rules 
for operators, with certain symmetry, sandwiched bet- 
ween symmetry-classified states. With our MSR forma- 
lism, one can of course recover the general selection rules 
predicted by WET, but in addition we can further ob- 
tain systematically simpler expressions for the amplitude 
of transitions through the help of UREF's. 

The WET gives selection rules for an operator A Va 
transforming like an irreducible representation T c bet- 
ween two states bearing the representation T; on the left 



and r r on the right : ( ip 



appear in the reduction of r c £gF r (or, in a symmetric way, 
if A\ does not appear in r ; * <g>r c <g>r r ). If A r " is a member 
of an ITO Ar£, the generalized Wigner-Eckart Theorem 
(gWET) gives more information about the amplitudes by 



factorizing ( rp 



tp^ r ) into the product of "reduced 



matrix elements" , which only depend on irreps Yi, r c , Y r 
and on physics-independent Clebsch-Gordan (CG) coef- 
ficients which involve also partner indices /ii,ii c ,fj, r [20| . 
With our formalism, it is easy to further separate the 

spinorial and spatial part of ( ip 



ijfc ) and sim- 



plify the result using only single group selection rules 
linked with UREF's. This procedure amounts to use fur- 
ther gWET selection rules at an intermediate level in the 
theory. 



A. Selection rules at the intermediate level of 
UREF's 

Equation (|39p is formally equivalent to write 
belongs to the OBB 



(r) 



Tb,ii},j where the ket 
The envelope functions 



^uf'SW °f sec tion HV] can be systematically decompo- 
sed into UREF's (/>~ i,rQ (r), which bear a single group 

irrep index r a and a corresponding partner function 
index fi a (c.f. Eq. (|igj)). 

Let us now return to the main problem of finding the 
matrix elements of ITO operators with, for simpli- 
city, the restriction that A^ operates only on spinorial 
coordinates. The matrix element can be written as 



€i\4:\<:) = 

^ a(Ti, m,Tb', (ib>;T r , /i r ,r&, 



x 



b' , A*f>' 



.4 



where the gWET can be 



used 
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to evaluate 



b' , fJ-b' 



r&, fJ-bj (this is nothing else than the 
matrix elements of a r c -ITO expressed in the OBB). 
The coefficients a(F;, m,Tv , fj,b']Y r , [i r ,Tb, Mfc) can be 
evaluated in turn as 



E 



a 



r t ,rr 



C, 



rv,r 6 * ; r Q 



6 r,r. 
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if does not denoted 



where we have used the gWET (at the intermediate level 

the so-called "reduced matrix elements" appear, they are 

fi,r a i ,rv,r a 



of UREF's) to evaluate 



As a result 



and defined by 



1 i. 



jr„r. 



l d r(cb T J> Ta (r 



(r) 

(93) 



Note that, according to the gWET, the reduced matrix 
elements are independent of the partner function index 
/j, a appearing in the right hand side. 

To conclude let us emphasize that such use of gWET 
and corresponding selection rules "at an intermediate 
level" is a salient feature enabled by the use of OBB and 
the resulting appearance of UREF's. Note that it allows 
in addition to minimize strictly the number of integrals 
to evaluate, furthermore the use of the SDR technique 
is enabled so as to compute all such integrals on the 
minimum domain. 

Besides these numerical advantages, we shall show in 
the example treated in the foregoing section that novel 
strong analytical features, physically observable, can also 
be most straighforwardly calculated with the use of the 
gWET "at an intermediate level" . 
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B. Polarization anisotropy for ground state 
momentum matrix elements in C^v symmetry 

Our example will be an operator important for optical 
transitions in semiconductors : the momentum operator 



P = P x e, 



P 



v °j» 



P z e z between a ground conduction 



band state and valence-band states in C$ v quantum wires 
and quantum dots (id (d = x, y, z) are unit vectors along 



the main polarization directions 22j, c.f. Fig. [7]). A small 
complication is that one must include here the conduc- 
tion band spin, by considering the products of the enve- 
lope function single group irreps A±,A 2 ,E with the irrep 
E1/2 for spin. The representation V_ B (g) corresponding 
to the OBB for valence band is given in Eq. (|42|) . For 
conduction band we choose V_ B (g) = D E (g)x Es/2 {g), 
corresponding to the central 2x2 block of V_ B (g). Since 
P x is Ai and (P V ,P Z ) form an ITO linked~with E, the 
application of the gWET at the intermediate level on 



b' , £ifc' 



Pi 



T b ,[ib),d = x,y,z leads to the following 



2x4 matrix representation for P x 



P P 



p. 


= Po 
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= Po 
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= Po 
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where P Q = (S\P X \X) = (S\P y \Y) = (S\P Z \Z) is the 
Kane matrix element. Another equivalent way to obtain 
them is to apply the U(c 2 ) transformation on the stan- 
dard Kane matrices 

Let us now consider the momentum operator matrix 
elements related to optical transitions from the valence 
band to ground conduction band state, assuming that 
the latter has A\ symmetry. This is generally the case, 
for common QWRs or QDs. Neglecting electron spin, we 
can use single group representations. Including electron 
spin, this ground state becomes twice degenerate and 



transforms like the product irrep A\ ® E- 



1/2 



E 1/2 , the 



conduction band spinor can then be written simply as 



— c,l 



(A, 



>(r) 



V c Al (r) 




^ 2{Al \r) 





^ (r) 
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where the index c recalls the nature of the state and may 
scan all conduction band states of symmetry A\. 

Using the conduction band spinors (|95t we shall eva- 
luate the squared momentum matrix elements between 
such E1/2 {A\) ground conduction band state and the 
v-th valence band state with irrep r„ (either 1 E 3 / 2 , 2 E 3 / 2 
or Ei/ 2 ). The squared momentum matrix elements are 
defined by 



M d (E l/2 (A x )-T v ) = Y, 



Pi 



ft: 



(96) 



where fj-c and fi v stand respectively for the partner func- 
tions of the conduction band state Ei/ 2 and the valence- 
band state r„. The scalar product can be decomposed 
as 



& (M) \P d 



/ , {E\/ 2 ,Hb>\Pd 



Ft, \ib 



x a(c : E 1/2 (Ai), He, E 1/2 , n b ,;v : T v ^ v ,T b , Hb) 

(97) 

where (-E1/2, fJ>b' \ Pd F b , fib*J and the coefficient a(c : 

Ei/2(Ai), He, E1/2, Hb>;v ■ T v ,Hv,F b ,Hb) are given by 
Eqs. (|94f and (|93|) respectively. 
The explicit results are as follows : 
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The first two results can be easily understood by recalling 
the fact that P x is an Ai-ITO : according to WET only 
r — r transitions are allowed since A\ ® T = Y. Then, 
coming back to a single group labelling for the conduction 



band, we find that the A\ — E- 



3/2: 



1, 2 transitions are 



forbidden in the a;-direction. 

However eqs. (|98|) contain a more striking result : from 
the second and fourth line we see that for the transition 
A\ —Exn there is a constant analytical ratio of oscillator 
strength (squared momentum matrix element) between 
the x- and y/z-directions. This entirely novel analytical 
prediction is in principle accessible to experimental verifi- 
cation, since the ground valence-band state (in the realis- 
tic QWR's investigated in [28[) is precisely oi E 1 / 2 sym- 
metry. It has the same value of polarization anisotropy as 
in quantum wells (22|, but it should be pointed out that 
this result is still a completely new result, much stronger 
than saying that the ground transition would be "quan- 
tum well light-hole-like" . Indeed this result holds exactly 
in presence of band-mixing as well, like at zone-center of 
the ground Ei/ 2 QWR subband, and also for any tran- 
sition with same initial and final symmetries, whether it 
would be in C 3v QWRs or C% v QDs, i.e. independently of 
dimensionality. In [48j a fairly detailed analysis of opti- 
cal polarization anisotropy in C% v QDs has been carried 
out, on the basis of the quantum well heavy /light-hole- 
like analogy, but the much stronger results obtained here 
are complementary and bring more light. The complete 
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analysis, including electron excited states as well as other 
kinds of hole states, will be reported elsewhere (28j . 

To summarize, we have illustrated in this section that, 
with the new formalism, the matrix elements of opera- 
tors do take a very simple analytical form, and that only 
overlaps of scalar UREF functions with the same symme- 
try are involved. Besides this new advantage, we do not 
expect additional selection rules at the global level. All 
the intermediate level calculations become minimal from 
a computational point of view, and much more trans- 
parent. It is the reason why we were able to find no- 
vel analytical ratio of oscillator strength between x- and 
y/z-polarization directions for certain optical transitions, 
which hold even in the presence of valence-band mixing. 



VIII. APPLICATION OF THE MSRF 
FORMALISM TO A VARIETY OF SYMETRIES 

Up to now we have discussed only the example of 
nanostructures with C% v symmetry. In this section we 
shall shortly discuss the application of MSRF forma- 
lism to other specific symmetries : C§ v , -D3/1, C n and 
C s . Cq v is a higher symmetry group presenting a 6 fold 
axis, some real QD structures like the wurtzite-based 
GaN/AIN QDs [3, |3tJ do display hexagonal symmetry 
like in Fig. [5] D^h symmetry can be found at zone-center 
in Cs v QWR's [28]. C n are pure rotational sub-groups 
(without roto-inversions), they typically correspond to 
the reduced symmetry occurring in previous structures 
under a magnetic field. Our last example will be the C s 
case, since it corresponds to the first type of structure - 
and the simplest - that we studied in the first section of 
this paper QTBJ] (see also [lCX |36|). 

In the following, we shall only discuss the symmetry of 
envelope functions obtained by taking into account the 
optimal choice of Bloch function basis and vectorial basis 
in real space. 



A. Ce, v symmetry 

One can explicitly construct Cq v from C$ v by adding 
a single 7r-rotation C 2 : C§ v — {g,gC2, Vg £ Cs v }. The 
single group C% v has four ID irreps Ai,Bi,i = 1,2, and 
two 2D irreps Ei , i = 1,2. Both the ID irreps and the 
two 2D irreps are half even and half odd with respect 
to the new operation Ci. For the double group, one has 
three 2D irreps E^n, i = 1,3,5. 

Let us now consider the p 7 q, r, s operators appearing 
in the valence band Luttinger Hamiltonian (c.f. subsec- 
tion |TTX|, and which are second order polynomials in k. 
We note P, Q, R, S the corresponding 3x3 "matrices" of 
coefficients in such way that p = k'Pk , ... , s = k'Sk. 
For C% v we recalled in subsection III Al that (r, s) form 
an ITO transforming with the E irrep, indeed it can be 
shown that the corresponding R and S matrices can be 



written as 

R o 3v 



a 
a b 
0-6 



S C 3v = 



0a 
-b 
a —b 
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where a and b are some spatially dependent constant de- 
pending on the Luttinger parameters linked with the un- 
derlying material. 

For Cq v further decompositions occur, and we can 
write r = ri + r 2 and s — s± + S2 where (ri,Si) are 
ITO's transforming like Ei irrep of Cq v . We obtain 



1 V 0/ 1 
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000 
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where a — Re (a) and b — ilm(b) are respectively associa- 
ted to the Ei and the E 2 ITO's. In Eqs. (TlTOf and ([TUT]) 
the effects of increasing the symmetry appear clearly : 1) 
the parameters a and b are simplified, 2) every parameter 
is related to a different irrep Ei in this new decomposi- 
tion. 

Another way to explicitly obtain the Cq v Luttinger Ha- 
miltonian starting from C% v would be to symmetrize the 
Hamiltonian with respect to the C 2 operation, as follows 

H c « (r, k) = 1 (£T C - (r, k) + Oa iT o »» (r, k)^ 1 ) 

(102) 

Of course we have to assume in addition that the spatially 
dependent Luttinger parameters in Eq. (|102[) are now 
invariant with respect to the Cq v symmetry operations. 

Switching to the OBB for Cq v , we find that the OBB 
reduces to Ex/2 ® -E3/2, an d that the Luttinger Hamilto- 
nian expressed in the 1/2, —1/2, 3/2, —3/2, basis (special 
order) reads 



H T = 



' (P-Q)h A 
m v A+ (P+Q)i2, 

(103) 

The decomposition of envelope functions into 
UREF's, expressed naturally with respect to the 
+3/2, +1/2, -1/2, -3/2 basis (standard order), leads 
to the following single group irreps (the argument r of 
functions is omitted in the following) : 
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(|T04|) indicate 



accor- 
U, that Ex/2 



The different components in Eqs. 
ding to the subduction tabie Cq v 

and E 5 / 2 are related to Ex/ 2 (C 3v ) while E 3 / 2 is related to 
1 E 3 / 2 ® 2 E 3 / 2 - As before reduced Hamiltonians can also 
be calculated. 



the corresponding irrep E i / 2 (D 3 h) , i — 1,5. In such a 
case (small symmetry breaking), pure Ds^-selection rules 
indicate which transitions have only very small matrix 
elements due to the approximate symmetry (see e.g. [28j| ) . 



D 3h symmetry 



The C n groups : subgroups of the rotations 
group 



The symmetry D 3 h is not common in semiconduc- 
tors nanostructures. However it can be shown in C 3v - 
symmetry Quantum Wires (QWR) that it is an approxi- 
mate symmetry at zone center [281 ] - It is obtained from 
C 3v by adding a symmetry operation ah, a planar re- 
flection with respect to an "horizontal" symmetry plane, 
i.e. perpendicular to the rotation axis of operations. 
Indeed D 3h = C 3v ® C s — {g,g<7 h , Vg € C 3v }, and the 
single group irreps of D 3h are simply the double of those 
of C 3v , and correspond to even and odd irreps with res- 
pect to the new horizontal symmetry plane Oh [2l| . For 
the double group irreps, the situation is sightly more com- 
plex : D 3 h involves then three 2D irreps Ei/ 2 , i = 1, 3, 5. 
The descent of symmetry tables [2l[ for D 3 y L — » C 3v give 
the correspondance E 3 / 2 -> 1 E 3 / 2 © 2 E 3 / 2 , E x /2 -> E\/2 
and E 5/2 Ex/2- 

In the case of the Luttinger Hamiltonian, the eigens- 
tate decompositions obtained with MSRF are as follows : 
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in the |j5 3 / 2 , l) , |£ 5 /2, l) , |#5/2, 2) , \E 3/2 , 2) basis. 

We shall naturally find additional selections rules in 
D 3 h 1 they are related to the different " D 3 h labels" cor- 
responding to otherwise similar Ex/2 eigenstates when 
C3i>-symmetry is only taken into account. Comparing the 
partner functions o f Ei/ 2 (D 3 h) and E 5 / 2 (D 3 h) presented 
in Eqs. ([106]) and (fT07|) to the Ex/ 2 (C 3v ) partner func- 
tions given in Eq. (|53j) . we note that the superpositions 
of Ai and E functions in the +1/2, —1/2 components, se- 
parate into two different irreps in D 3 h (only even scalar 
function respect to <j% are involved). 
When the deviation from D 3 h symmetry is very small, 
every Ex/ 2 {C 3v ) gets only a small A\ or E part in agree- 
ment with the partner function decomposition linked to 



The subgroups C„ , n £ N of the pure rotation group 
SO (3) result from a descent of symmetry from C nv . As an 
illustration we compare in Fig. [T5] (a) and (b) the shape 
of a C/x v and a C4 quantum structure. They are also of 
particular interest in the important case of an applied 
magnetic field on a C nv structure (see e.g. the case of 
In As quantum dots [43]). 



(a) 



(b) 



Figure 13: General shape of quantum structure with sym- 
metry, (a) C iv ; (b) C 4 

Since C n groups are cyclic and abelian groups, all their 
single and double group irreps are one-dimensionnal. 
Therefore in the C n case, the reduction of representations 
is equivalent to a diagonalization of the representations. 
This makes rather trivial the approach in (4^ . [50| , where 
the authors study C^ v and Cq v QDs with and without 
magnetic field using only the properties due to C n sym- 
metry. To compare with the present work no additionnal 
separation in spatial and spinorial representation is per- 
formed, there is only a block-diagonalization of the full 
Hamiltonian into ID irreps of C n before numerical com- 
putation. We shall show that even in such a simple C n 
case, the use of a fully symmetrized basis together with 
spatial reduction sheds some new light and allows some 
further analytical and numerical simplifications. 

Let us study the C4 symmetry, which is schematically 
represented in Fig. I13l bv comparison with C^ v . The four 
single group irreps of the group C4 are A,B, l E, whilst 
its four double group irreps are i E 1 /2, l E 3 / 2 , i = 1, 2. The 
l E, l Ex/ 2 , % E 3 / 2 irreps with % = 1,2 form three sets of mu- 
tually conjugated one-dimensional irreps. Therefore in 
the OBB the V_ B (g) representations matrices are neces- 
sarily diagonal and would correspond to 1 £ , 3 / 2 © 1 -Ei/2 © 
2 Ex/ 2 © 2 E 3 / 2 in the case of our Luttinger Hamiltonian. 
By contrast, one can see that, despite the fact that the 
3x3 representation matrices constructed with the ro- 
tation matrices (|12p are block-diagonal, these blocks are 
not naturally labelled by irreps of C4. The lxl block is 
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trivially linked with the irrep A, and the reduction of the 
2x2 block leads to 2 E © 1 E. The required switch from 
the cartesian vectorial basis to the "Optimal Vectorial 
Basis" (OVB) is completely analogous to the OBB, and 
is given by 



e A = e. 



1 

i 



(e„ - ie z ) 



(108) 



The last two vectors correspond, in optics, to complex 
unit vectors for circular polarization in the y — z plane. 
The OVB has been used in ref. [45j to treat a vectorial 
problem linked with the eigenmodes of photonic crystal 
microcavities. 

Using such fully symmetrized OVB, simultaneously 
with the OBB basis, one can also obtain for C4 the 
ultimate decomposition of spinorial eigenstates into 
UREF's ; 



^,£3/2 



(109) 



As we now would like to discuss selection rules for opti- 
cal transitions between the conduction and valence band, 
one should remark that the use of the "optimal vectorial 
basis" calls for new corresponding Kane matrices. To this 
end a preliminary remark is that for conduction band 
states with spin (j = 1/2), the V_ B (g) representation re- 
duces to diagonal matrices according to 1 E 1 / 2 © 2 ^i/2- 
We also note that, in a similar way as for C% v , the spi- 
nors decomposition correspond to the central part of the 
corresponding valence band decomposition (j = 3/2). 

The new Kane matrices are then 
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As far as selection rules are concerned, it can be 
checked term by term with Eqs. (jl09[) and (|110[) that 
the allowed transitions between a conduction and a va- 
lence band state follow indeed the C4 irrep multiplica- 
tion table [2TI ] . i.e. only diagonal r — T transitions are 
permitted in x direction (where T is any irrep), and 
l E n - (3 ~ l) £ (2 -n) (with i = 1,2 and n = 1/2,3/2) are 
only allowed in er_ polarisation, whilst only E n — ( 3 ~ l 'i?„ 
are only allowed in er + polarization. 



D. Return to the C a group 



Although we developped the new formalism to study 
the more difficult HSH, it is interesting to return to 
C s low symmetry heterostructures, like T or V-shaped 
QWRs seen in III B II and treated by the Luttinger Ha- 
miltonian, to see how the novel formalism reduces to the 
old solution technique : the choice of an " optimal quanti- 
zation axis" perpendicular to the symmetry plane a [l(| 
(see also [Hj]). This will also definitely put all cases on 
the same firm ground. 

The "optimal quantization axis choice" was only a 
way to choose the OBB which diagonalized the V_(a) 
related to the planar reflection with respect to the a 
symmetry plane. With this we obtained even/odd en- 
velope functions [To| which we can now associate to 
UREF's labelled by the one-dimensional single group ir- 
reps A', A" . However we see that the complete reduction 
of the 4x4 matrix representation V_(g) of C s may read 

as 1 E 1 / 2 © 2 E 1 / 2 © 1 £7 1 / 2 © 2 ^i/2- A few remarks are here 
in order : 1) since all double group irreps of C s are one- 
dimensional it amounts, this particular case, to a diago- 
nalization of the representation (like in C4), 2) the uni- 
tary transformation is accomplished by a rotation matrix 
parametrized by the three Euler angles (a, /?, 7) corres- 
ponding to a "3D rotation" 3?(a,/3,7) of the so-called 
quantization axis direction, 3) j — 3/2 labels can be kept 
for the new basis and refer to rotated Bloch functions, 
4) we see that for the first time some irreps appear twice 
in the decomposition ( 1 E 1 / 2 and 2 E 1 / 2 ). Considering the 
cartesian vectorial basis presented in Fig. [71 and a verti- 
cal plane normal to the e z direction, the vectorial basis 
e x ,e y , e z is already optimal and reduces to A' © A' © A" . 
To summarize, the Luttinger Hamiltonian eigenstates can 
be decomposed in the basis above as 



n/2 



V 



(111) 



where, for every ^1/2 irrep, (f) V and (p T are different func- 

1— 

tions with the same symmetry. The functions in ip 



'1/2 



and ip 1/2 are identical only if ijj 1/2 and ip 1/2 are rela- 
ted by time reversal symmetry. As already seen for C4, for 
the conduction band with spin, the j = 1/2, m = ±1/2 
envelope functions would just be analogous to the cen- 
tral part of the j = 3/2, m = ±1/2 valence band enve- 
lope functions. The momentum operators P x and P y in 
the symmetry plane are even with respect to a (A') and 
the operator P z perpendicular to the axis is odd (A"). 
The latter A" operator only couples mutually conjuga- 
ted bands (M^' r = My :F = 0, selection rule) whilst the 
former A' operators only allow diagonal T — T transitions. 
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IX. WIDE POTENTIAL APPLICABILITY AND 
ORIGINALITY OF THE MSRF FORMALISM 

The foundation of our new MSRF formalism can be 
summarized, for most general spin dependent problems, 
as first choosing a really optimal basis (fully symmetri- 
zed basis) for the orbital and spinorial space and trea- 
ting every spinorial component as a sum of symmetrized 
UREFs. Second, on this set of UREFs the spatial do- 
main reduction (SDR) technique can be applied. With 
the fully symmetrized basis, the coupling between dif- 
ferent envelope functions becomes minimized and every 
spinorial function can be decomposed with respect to 
symmetry using only single group irreps. With the SDR 
applied to every UREFs one can then achieve domain re- 
duction. In this context it should therefore be clear for 
the reader that in band-structure problems any number 
of bands could easily be treated by the MSRF at the 
price of more complexity, strain [4JJ could also be trea- 
ted (treating also by MSRF the separate strain equations 
in a self-consistent way) , as well as the Burt-Foreman in- 
terface terms [23|, [5l|, [13] , or the presence of an external 
field like a magnetic field (if one takes into account the 
lower global symmetry in this case). The choice of fully 
symmetrized basis presented in Sec. HVI onlv depends on 
the symmetry group considered (some other examples 
have been given in Sec. I Villi) . 

The SDR technique used in Sec. |V] can be easily ge- 
neralized, but one must realize that it may need slight 
tuning depending on the problem at hand. For example 
for higher order FEs (like second order FEs which are 
frequently used), the general SDR procedure would be 
exactly the same, provided one is careful in considering 
anew the connectivity between sub-domains (a few new 
blocks must be included). 

It is important to note that the MSRF formalism is not 
restricted to our example of a 2D k.p spinorial problem 
but can be easily generalized for any other more general 
vectorial, spinorial or tensorial-like problems, also with 
different dimensionality (1D-3D), and a numerical solu- 
tion with other techniques, whether they would be in the 
spatial domain or the Fourier domain, or in any other sen- 
sible basis. The approach is even not limited to linear pro- 
blems provided one accounts for proper products of repre- 
sentations. The essence is to symmetrize separately the 
bases in the spinorial-like space and the real "coordinate- 
like" space. Real space resolution methods, are usually 
the most powerful for the study of spatially isolated he- 
terostructures (i.e. occurrence of band matrices or sparse 
matrices), whilst Fourier space decomposition of the en- 
velope functions, is most powerful for periodic structures 
(i.e. superlattices). Indeed for a non-periodic problem a 
formulation in a plane wave basis, which give rise to full 
matrices during the numerical solution, would lead to a 
significant waste of computer time and memory space 
compared to real space methods like FD of FE methods. 

In real space formulations the MSRF method sketched 
above has the further potential to keep the sparse ma- 



trix structure of the numerical problem (when formula- 
ted with a FD or FE approach) , therefore it can be mar- 
ried with convenient and maximally efficient methods for 
sparse matrices. 

The MSRF formalism is also plainly applicable to the 
study of strongly coupled periodic structures where sym- 
metrized plane waves decomposition in the way of Vuk- 
mirovic et al. [49[ are efficient. For a k-p Hamiltonian 
the optimal Bloch function basis would of course be the 
same, the only part that has to be modified is the spa- 
tial domain reduction approach, where one would then 
introduce the concept of a " reduced plane waves basis" . 
We have already achieved steps in this direction in [44^ . 
where a procedure to work with a reduced domain in Fou- 
rier space has been introduced (see Fig. 6c of 44] to sec 
the reduced Fourier domain for every irreps of the C^ v 
group). Once this is done the way to proceed for stron- 
gly coupled periodic structures is along the same lines, 
and would again lead to optimal algorithms and detailed 
symmetry analyses. 

It is also worth to point out that once a calculation is 
done the completely symmetrized envelope functions pro- 
vided by the MSRF formalism may also be very useful 
to further build in a maximally efficient way more com- 
plex symmetrized objects. A typical example in semicon- 
ductor heterostructures in building excitons, or excitonic 
complexes, which gives rise to multidimensional Coulomb 
integrals involving products where not only additional 
single group selection rules can be used, but where also 
the domain reduction would allow drastic speedups of 
the computation of integrals. Of course in this domain 
too MRSF also has the potential to provide further na- 
tural physical insight by the resulting minimal coupling 
scenario at the individual envelope function level. 

To which extent does the MSRF links with previous 
work using heterostructure symmetry to simplify the cal- 
culation of electronic levels ? Little work has in fact been 
done up to now on this topic. Let us compare our method 
to the work of Vukmirovic et al. 143. 15CI which appears 
to be closest of ours, in particular in [49(, a method using 
symmetrized plane waves decomposition is suggested for 
a Ci pyramidal QD. Our first comment is that such a 
Fourier decomposition is most useful only for strongly 
coupled periodic structures, like found for example in [50| 
for hexagonal QD superlattices. In the other limit of spa- 
tially isolated heterostructures the approach is not nume- 
rically optimal because it leads to full matrices instead 
of band matrices like in real space approaches like FD or 
FE schemes. Our second comment is that symmetry is far 
from fully exploited in their approach. Let us discuss this 
last point in more details. First the authors define an N- 
dimensional representation of the rotation operators on 
the full coupled Fourier-spin space (no separate symmc- 
trization), and they show that a "symmetry- adapted ba- 
sis" can be found (diagonalization of the representation) . 
As expected this allows to block-diagonalize the Hamilto- 
nian, but really it should be pointed out that the block- 
diagonalization achieved is nothing else than the main 
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decoupling that one should obtain by separating the ei- 
genproblem into its main different solutions, indexed by 
double-group irreps. Second, the authors show that their 
method allow to separate the Hamiltonian into n blocks 
of equivalent sizes which does reduce the necessary com- 
putation time, but this is obvious for a C n group. The 
limitations of the approach [49|,[5(| are therefore evident. 
First, it does not separate orbital and spinorial contribu- 
tions (only double group irreps appear in the reduction) , 
hence the symmetry properties of the structure are ta- 
ken into account only at the last step before numerical 
resolution, no single group classification of envelope func- 
tions, no selection rules at the envelope function level can 
be obtained, nor any domain reduction can be achieved. 
Second, this approach is well-adapted only for the case 
of rotation sub-groups C n where every irrep of the single 
and double group are ID. As we have seen earlier already 
for small dihedral groups like Civ > degenerate eigenvalues 
could appear. Indeed in |50j the author study a Cq v pro- 
blem (every double group irreps is 2D) but only exploit 
the Cq symmetry (with only ID irreps) by neglecting the 
vertical symmetry planes, since in their framework they 
cannot describe optimally the further degeneracies of the 
eigenstates. Of course in the presence of a magnetic field 
this would be the correct approach, since additional lif- 
ting of degeneracies would occur. 



X. CONCLUSION 

We have presented a new formalism called MSRF 
(Maximal Symmetrization and Reduction of Fields) wi- 
dely applicable to the study of high symmetry hete- 
rostructures of various dimensionality (e.g.QWR's or 
QD's), and well adapted to treat both scalar and spinorial 
problems. The use of this formalism has been illustrated 
with particular examples, mainly a C^v QWR. 

For scalar functions (e.g. conduction band in the frame 
of a single band spinless Hamiltonian) it reduces to the 
spatial domain reduction (SDR) technique. For every ir- 
reducible representation (irrep), the independent para- 
meters provide systematically the reduced minimal do- 
main for significant description of every quantum state, 
and allows to throw away the computation of all the re- 
dundant parts by identifying the reduced Hamiltonian 
for the independent parameters. The method was deve- 
loped in a pedestrian way for C$ v , further mathematical 
developments of the SDR technique will be reported el- 
sewhere. 

For spinorial sets of functions (e.g. valence band with 
a four bands Luttinger Hamiltonian), it is necessary to 
separate the spinorial and orbital parts and to compute 
a spatial domain reduction of the spinorial problem. For 
this purpose we introduced the concept of optimal fully 
symmetrized basis for the spinorial field (every basis 
member transforms like a partner function of an irrep 
of the group). The first advantage is the simplification 
in spin space, since the Optimal Bloch function Basis 



(OBB) completely reduces the matrix representation of 
the transformation laws in spin space into block-diagonal 
form, which minimizes the complex coupling between spi- 
norial components by symmetry operations. Moreover we 
showed that every spinorial component can be decompo- 
sed further into ultimately reduced envelope functions 
(UREFs), labelled by single group irreps, which finally 
allows the application of the SDR technique. 

The advantages of such a formalism are manifold. In 
particular the Hamiltonians take usually a much more 
simple form. Operators and spinorial component can sim- 
ply be labelled with single group irreps, and with the 
help of the SDR technique one can solve a smaller speci- 
fic problem for every irrep on a reduced spatial domain. 
This procedure leads to simpler analytical expression for 
operator matrix elements. Just as with any symmetrized 
formalism, any symmetry breaking mechanism can be un- 
derstood more easily, both in a qualitative and quanti- 
tative way. Hidden approximate higher symmetries can 
also be quickly detected by the analysis of selection rules. 
From the numerical point of view the formalism leads to 
highly advantageous algorithms, and applies in a flexible 
way to a wide variety of methods (real or Fourier space), 
so that the most practical one for the case at hand can 
be chosen without restriction. The reduction factor in 
computer time for diagonalization can be estimated to be 
about 20 for C^ v . Last, but not least the MSRF technique 
can be applied as a post-symmetrization technique on 
existing numerical results (no absolute need to recode), 
giving access to the full power of symmetry analysis at 
the envelope function level, and greatly increased perfor- 
mance in subsequent numerical computations. In |45l.l47j 
we have implemented the automatic recognition of mode 
symmetries, as well as the postsymmetrization technique. 

In a forthcoming paper [28j extensive analytical and 
numerical results on the electronic and optical properties 
of a real C% v QWR will be presented and will further 
demonstrate the power of the MSR approach. 

As already stressed, straightforward generalizations of 
the method may be developed for arbitrary tensorial 
fields obeying a set partial differential equations (even 
non-linear!). Therefore in a forthcoming paper [45[ we 
shall apply it to Maxwell's equations in a case corres- 
ponding to photonic bandgap microcavities [53|. Other 
possible application could be connected heterostructure 
problems (definition of strain fields, phonon fields, etc.). 

A further high potential of MSRF also lies in possible 
subsequent calculations. Often, in a second stage, one is 
also interested to build more complex objects, like exci- 
tons, polaritons, polarons, etc... To build such objects the 
symmetrized field components provided by our technique 
really represent optimal building blocks, for which selec- 
tion rules and well-defined transformation properties are 
readily available by construction, hence MSRF should be 
of great use to an even larger community. Last, but not 
least, it is worth pointing out that the spirit of the MSR 
approach can also be applied to other widely used models 
for heterostructures like tight-binding, pseudo-potential, 
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etc., with the same potential analytical and computa- 
tional advantages. 
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Annexe A: C 3v point group tables 

In the first part of this appendix we recall the 
most important C$ v point-group tables which are 
extensively used in the main body of the paper (ad- 
ditional informations and tables can be found in |2l| . 
note that the multiplication table (c) is transposed with 
respect to ref . [2l[ since we use the passive point of view) . 
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Table I: Character (a), direct product of irreps (b), multi- 
plication (c) and factor (d) tables for C 3v group 



In the second part, since infinitely many equivalent ma- 
trix representations can be used for the 2D degenerate 
irreps of C3,,, we explicitly specify which 2D matrix re- 
presentations are used. For the E irrep we have chosen 
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the following 2D matrix representation 
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For the other 2D £1/2 irrep, the chosen 2D 
sentation can be constructed with the help 
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matrix repre- 
of the choice 



[1] C. Girard, Rep. Prog. Phys. 68, 1883 (2005). 

[2] A. D. Yoffe, Advances in Physics 51, 799 (2002). 

[3] K. Sakoda, Optical properties of photonic crystals (Sprin- 
ger, Berlin, 2005). 

[4] R. M. Stevenson, R. J. Young, P. Atkinson, K. Cooper, 
D. A. Ritchie, and A. J. Shields, Nature 439, 179 (2006). 

[5] R. J. Young, R. M. Stevenson, A. J. Hudson, C. A. Nicoll, 
D. A. Ritchie, and A. J. Shields, Phys. Rev. Lett. 102(3), 
030406 (2009). 

[6] S. Kako, K. Hoshino, S. Iwamoto, S. Ishida, and Y. Ara- 

kawa, Appl. Phys. Lett. 85, 64 (2004). 
[7] A. D. Andreev and E. P. OReilly, Phys. Rev. B 62, 15851 

(2000). 

[8] H. Akiyama, J. Phys. : Condens. Matter 10, 3095 (1998). 
[9] E. Runge, Solid State Phys. : Advances in Research and 
Appl. 57, 149 (2002). 
[10] M.-A. Dupertuis, E. Martinet, D. Y. Oberli, and E. Ka- 

pon, Europhys. Lett. 52, 420 (2000). 
[11] M.-A. Dupertuis, Phys. stat. sol. (b). 221, 323 (2000). 
[12] F. Vouilloz, D. Y. Oberli, M.-A. Dupertuis, A. Gustaf- 



sson, F. Reinhardt, and E. Kapon, Phys. Rev. B 57, 
12378 (1998). 

[13] P. C. Sercel and K. J. Vahala, Phys. Rev. B 44, 5681 
(1991). 

[14] P. C. Sercel and K. J. Vahala, Appl. Phys. Lett. 57, 545 

(1990) . 

[15] U. Bockelmann and G. Bastard, Europhys. Lett. 15, 215 

(1991) . 

[16] D. S. Citrin and Y. C. Chang, Phys. Rev. B 43, 11703 
(1991). 

[17] A. Hartmann, Y. Ducommun, K. Leifer, and E. Kapon, 
J. Phys. : Condens. Matter 11, 5901 (1999). 

[18] F. Michelini, M.-A. Dupertuis, and E. Kapon, Appl. 
Phys. Lett. 84, 4086 (2004). 

[19] P. Tronc, V. P. Smirnov, and K. S. Zhuravlev, phys. stat. 
sol. (b) 241, 2938 (2004). 

[20] J. F. Cornwell, Group Theory In Physics, Vol I (Acade- 
mic Press, London, 1984). 

[21] S. L. Altmann and P. Herzig, Point-Goup Theory Tables 
(Clarendon Press, Oxford, 1994). 



33 



[22] G. Bastard, Wave Mechanics Applied to Semiconductor [40 

Hetero structures (Les Editions de Physique, Paris, 1988). 
[23] B. A. Foreman, Phys. Rev. B 48, 4964 (1993). 
[24] B. Lassen, V. C. L. Y. Voon, M. Willatzen, and R. Mel- [41 

nik, Solid State Comm. 132, 141 (2004). 
[25] O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 

59, 5688 (1999). [42 
[26] P. Boucaud and S. Sauvage, C. R. Physique 4, 1133 [43 

(2003). 

[27] Q. Zhu, E. Pelucchi, S. Dalessi, K. Leifer, M.-A. Duper- 

tuis, and E. Kapon, Nano Letters 6(5), 1036 (2006). 
[28] S. Dalessi, F. Michelini, and M.-A. Dupertuis, Phys. Rev. [44 

B (2009), in preparation. 
[29] J. M. Luttinger and W. Kohn, Phys. Rev. 97, 869 (1955). [45 
[30] J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 
[31] G. Fishman, Phys. Rev. B 52, 11132 (1995). [46 
[32] A. Messiah, Mecanique quantique (Dunod, Paris, 1969). 
[33] F. Bassani and G. P. Parravicini, Electronic states and 

optical transitions in solids (Pergamon Press, Oxford, [47 

1975). 

[34] M.-A. Dupertuis, D. Marti, and F. Michelini, phys. stat. 

sol. (b) 234, 329 (2002), note : in this paper there are [48 
two missprints : in Eq. (1) W and W~ x are inverted and 
the same with K in Eq. (10). 

[35] D. H. Marti, M.-A. Dupertuis, and B. Deveaud, Phys. [49 
Rev. B 72, 075357 (2005). 

[36] D. H. Marti, M.-A. Dupertuis, and B. Deveaud, IEEE J. [50 
Quant. Electron 41, 848 (2005). 

[37] D. Simeonov, E. Feltin, J.-F. Carlin, R. Butte, M. lie- [51 
gems, and N. Grandjean, J. Appl. Phys. 99, 83509 (2006). [52 

[38] W. Appel, Mathematiques pour la physique et les physi- [53 
ciens (H & K, Paris, 2002). 

[39] M. Carmeli and S. Malin, Theory of spinors : an intro- 
duction (World Scientific, Singapore, 2000). 



E. P. Wigner, Group theory and its application to the 
quantum mechanics of atomic spectra (Academic Press, 
New York, 1971). 

G. L. Bir and G. E. Pikus, Symmetry and Strain-induced 
Effects in Semiconductors (John Wiley and Sons, New- 
York, 1974). 

A. Bossavit, SIAM J. Appl. Math. 53, 1352 (1993). 

A. Faessler and E. Stiefel, Group theoretical methods 
and their applications (Birkhauser Boston, Boston MA, 
1992), with English translation by Baoswan Dzung 
Wong. 

D. Obreschkow, F. Michelini, S. Dalessi, E. Kapon, and 
M.-A. Dupertuis, Phys. Rev. B 76, 035329 (2007). 

B. Gallinet, M. A. Dupertuis, J. Kupec, and B. Witzig- 
mann (2009), in preparation. 

R. D. Cook, D. S. Malkus, M. E. Plesha, and R. J. Witt, 
Concepts and Applications of Finite Element Analysis 
(John Wiley And Sons Ltd, United States, 2001). 

D. Dupertuis, B. Gallinet, and M. A. Dupertuis 
(2009), A postsymmetrization code will be posted on : 
|http : / /wiki.ep fl.ch/msrf-sdr[ 

K. F. Karlsson, V. Troncale, D. Y. Oberli, A. Malko, 

E. Pelucchi, A. Rudra, and E. Kapon, Appl. Phys. Lett. 
89, 251113 (2006). 

N. Vukmirovic, D. Indjin, V. D. Jovanovic, Z. Ikonic, and 

P. Harrison, Phys. Rev. B 72, 75356 (2005). 

N. Vukmirovic, Z. Ikonic, D. Indjin, and P. Harrison, J. 

Phys. : Condens. Matter 18, 6249 (2006). 

M. G. Burt, J. of Phys. : Cond. Mat. 4, 6651 (1992). 

M. G. Burt, J. of Phys. : Cond. Mat. 11, R53 (1999). 

O. Painter and K. Srinivasan, Phys. Rev. B 68, 035110 

(2003). 



